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ABSTRACT 


An expanded package of interactive FORTRAN computer pro- 
grams has been developed for the analysis and forecasting 
of time series data. The package, called the Time Series 
Editor, is designed to employ the iterative Box-Jenkins 
methodology of time Series analysis. The Time Series Editor 
was developed for time-shared use on the Control Program/ 
Cambridge Monitor System (CP/CMS) at the U.S. Naval Post- 
graduate School, but can be modified for use on other time- 
sharing systems with a FORTRAN capability. The Time Series 
Editor assists in data preparation and entry, analysis, 
modeling, forecasting and diagnostic testing. Utilization 
of the package, following the included User's Guide, requires 
only a limited knowledge of the computer system, with all 


required user responses interactively prompted by the Editor. 
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I. INTRODUCTION 


Operations researchers, military planners and programmers, 
statisticians, es marketing personnel, managers and 
Pavoeoeace OLten Laced with the meed Co analyze data in the 
form of a time series, which can be thought of as a Sequence 
Pemowservatilons Of either a deterministic or stochastic 
Biecess. In most cases, the objective of the analysis is 
to determine patterns or other recognizable behavior appar- 
ent in the data, and to then formulate a suitable mathemati- 
cal model for the time Series from which forecasts of future 
behavior can be obtained. The value to a decision maker 
of being able to predict the future with some reasonable 
and statistically quantifiable degree of accuracy can not 
be overstated. For example, in areas such as budget expen- 
ditures, recruiting performance, commodity prices, population 
levels, resource consumption, manpower levels and consumer 
demand for products, decision makers charged with planning 
for the future should base their decisions at least in 
part on the best available predictions about the future 
behavior of the time series in question. 

Until the late 1960's, the techniques employed in time 
series analysis were primarily those of Spectral analysis, 
with applications of harmonic analysis and mathematical 
transform theory. For an introduction to these topics, the 
reader is referred to Anderson [Ref. 2], Lewis [Ref. 7] 


and Jenkins [Ref. 10]. Due to the high degree of mathematical 





sophistication required by the spectral analysis approach, 
the capability to perform time series analysis resided nearly 
exclusively with mathematigians and electrical engineers. 
As a consequence, the majority of decision makers came to 
use more understandable (and far less powerful) methods 
such as simple moving averages or exponential smoothing. 
However, in the late 1960's and early 1970's, the statistical 
analysis of time series by the methods developed by Box 
and Jenkins [Ref. 4] has gained widespread acceptance. 
Although these methods are mathematically nearly equivalent 
to the spectral approach (there are transformations that 
interconnect the two methodologies), the Box-Jenkins approach 
is described in vocabulary more familiar to operations 
researchers, statisticians, economists and managers, and 
is therefore being used more and more by them in building 
models from measurements of the past. 

Many algorithms and computer programs for performing 
the analyses required by the Box-Jenkins approach have been 
developed and are available from several sources. Perhaps 
the best source is the collection of FORTRAN computer sub- 
routines which resides in the International Mathematical 
and Statistical Library (IMSL) [Ref. 9]. The primary prob- 
lem with using the available computer resources lies not 
in any deficiency of the programs themselves, but with the 
very nature of the Box-Jenkins methodology itself. The 


Box-Jenkins method is an iterative approach to time series 
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analysis, which is briefly described in Figure l. [See 
Wheelwright and Makridakis, Ref. 16.] 

As indicated by Figure l, the Box-Jenkins method is a 
multi-stage, iterative process. It begins with the postu- 
lation of a general class of time series models which has 
been found, experimentally, to be extremely rich. There- 
after, the procedure continues as a trial-and-error process, 
with decisions points where the analyst is required to 
select the next direction based on the best information 
available to him. Since each stage of the process outlined 
in Figure 1 may consist of several sub-steps, the modeling 
process itself can become quite time-consuming, even with 
the ready availability of the IMSL software resources. For 
example, an analyst employing the IMSL subroutines in a 
batch processing computer system to perform Box-Jenkins 
modeling of a time series might perform the following 
sequence of tasks: 

1. Prepare the time series data in the proper format. 

2. Plot and visually examine the time series, checking 
for nonstationarity, trends, seasonality, patterns, 
etc. 

3. Write a program to call the IMSL subroutine that 
will calcuate the mean, variance, autocorrelations 
and partial autocorrelations of the series. 

4. Plot the autocorrelations and partial autocorrela- 
tions; this provides much of the information re- 
quired for identifying the correct model for the 


series. 


LAL 










SELECT ANOTHER 
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ssa SGinew | 
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5. Write a program to call the IMSL subroutine which 
transforms the time series to adjust for seasonal 
patterns, nonstationary behavior or other behavior 
which deviates from that assumed for the class of 
models postulated. 

6. Repeat steps (2) through (4) with the transformed 
Glave. 

7. Review the statistical properties of the auto- 
correlations and partial autocorrelations for 
tentative model identification. 

8. Write a program to call the IMSL program that 
estimates the model parameters and computes the 
residuals. 

9. Write a program to perform goodness-of-fit tests 
for the model. 

10. Analyze the model residuals, using steps (1) through 
(9), as with the original time series. 

ll. Refine the model using information gained by 
examination for any structure in the residuals. 

12. Repeat the preceding process until no structure 
remains in the residuals. 

13. When an adequate model has been obtained, write a 
program to call the IMSL subroutine which forecasts 
and determines confidence intervals for future 
values of the time series. 

Between each succesSive pair of steps, the analyst must 


manually intervene in the analysis process, making a decision 


Le 





subjectively based on the information then available. 

Thus, much analyst interaction iS necessary in order to 
determine a suitable mathematical model and forecast equa- 
tion. Even with rapid computer job turnaround time, the 
process described above could easily consume a full working 
day or more. 

This paper describes an effort to anne aes some of the 
organizational problems inherent in ne node ling of time 
series using the Box-Jenkins iterative approach. This 
effort consists of an interactive computer program package 
which provides access in a structured way to the most useful 
IMSL and other subroutines necessary for Box-Jenkins time 
series analysis. This package, called the Time Series Editor, 
1S written for use with the Naval Postgraduate School's 
Control Progam/Cambridge Monitor System (CP/CMS). Since all 
working programs with the exception of the executive routine 
itself are written in FORTRAN, the Time Series Editor should 
be adaptable for use on other FORTRAN-capable time-sharing 
systems. The Time Series Editor assists the analyst in 
data preparation and entry, model construction and diag- 
nostic testing, and time series forecasting. In fact, 
with the User's Guide provided as an Appendix to this report, 
a complete Box-Jenkins time series analysis can be per- 
formed in a short time by even a novice computer user. 

For this reason, the Time Series Editor could be a valuable 
instructional tool for laboratory use in a time series 


eOurse. 
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A brief overivew of Box-Jenkins methodology is provided 
in Chapter II, to serve as a point of reference for the 
material which follows. Chapter III contains a description 
of the algorithm employed to estimate non-linear least- 
squares parameters for a generalized (seasonal, nonseasonal, 
Stationary, or nonstationary) Box-Jenkins model. Chapter IV 
contains descriptions of each of the major programs in 
the Time Series Editor. Chapters V and VI contain examples 
of actual use of the Time Series Editor, using a non- 
seasonal and a seasonal series, respectively. Chapter VII 
includes a summary of the report, and recommendations for 
future additions to the Time Series Editor. The Appendices 
include a User's Guide to the Time Series Editor, with 
instructions, sample user sessions and sample outputs, as 
well as listings of the Editor itself and all primary 


Subprogram packages. 


nes: 





II. BOX-JENKINS METHODOLOGY: A PRIMER 


This chapter provides a brief overview of the methodology 
developed for the analysis of time series data by G.E.P. 
Box and G.M. Jenkins, in order to facilitate understanding 
of the models and programs employed in the Time Series 
Editor. For more detail on the material covered in this 
Saapter, the interested reader is referred to the texts by 
Anderson [Ref. 1], Box and Jenkins [Ref. 4], Pindyck and 


Rubinfeld [Ref. 15], Nelson [Ref. 14], and Mabert [Ref. 12]. 


A. LISTING OF NOTATION 
The following is a listing of the notation employed 
in this chapter and in the remainder of this report. Each 


symbol is provided with a brief title or explanation of its 


use. 
{z,} a time series 
N the length of the time series (number of terms) 
Zi, 5? _ estimates of the mean and variance of the time 
5 series {z,} 
r 
’ ee = 
{zi} transform ene {zi}, where Zo (2, e) 
' = + 
or Zz) In(z, 5) 
yt differenced series from {fz}, where 
Ww = yous. 
ce s't 
Py autocorrelation of lag k 
ry estimate of autocorrelation at lag k 
d 
V backward difference operator of order d 
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seasonal backward difference operator of order 
D with season s 


backward shift operator of order s 
length of the time series' season 


non-seasonal autoregressive parameters 
(estimates, $) 


seasonal autoregressive parameters (estimates, 6) 


non-seasonal moving average parameters 
(estimates, 6) 


seasonal moving average parameters (estimates, 0) 
series of white noise terms of the process 
variance of the white noise process 

i 

sum of squared residuals of the model 

number of non-seasonal autoregressive parameters 
number of seasonal autoregressive parameters 
number of non-seasonal moving average parameters 
number of seasonal moving average parameters 
transfer function of a linear filter 


partial autocorrelation of log k 


ae ME SERIES, STOCHASTIC MODELS AND STATIONARITY 


A distinguishing feature of modern time series analysis 


is that the sequence of observations of a given variable is 


considered to be a realization of jointly distributed 


ai 





random variables, such that the time series itself can be 
viewed aS a stochastic process. A discrete stochastic 
yizgre ey) 


generated sequentially in time by a set of jointly distributed 


time series, then, 1S a set of observations (Zz 


random variables, where the actual data (Zp 1Zor--- 72) 


represents a particular realization of some joint probability 


wrstribution function £(z Wie Ches funet Lon 


a 


determined in some way, a forecast Zi 4, Can be thought of as 


having been generated by a conditional probability distri- 


bmeron f(z ). Now, the stochastic process 


pel 2prZgre rey 


which generates the time series (Zy1Zore-- 2 ) 1S said to 


is 


be stationary if its properties are unaffected by a change 

in time origin; that is, the process has reached a particular 
State of statistical equilibrium. This implies that the 
joint probability distribution associated with the m 


observations (Zz ), Made at any set of times 


ea eaey A 
al 5 th 


(ty,t5,---,t), is identical to that distribution associated 


with the m observations (2z ) made at 


PayA peeez 
t,+k to+k t tk 


times (t,+k,t +k,..-+,t +k). 


1 2 
For stationary discrete time series, the marginal 
Beobability distribution p(z,) is the same for all times 


t. Hence, the stochastic process embodied in the time 


series has a constant mean 


y= Elz] = I SZeeay azo; 


ts: 





which defines the level about which it fluctuates in time, 


and a constant variance 
2 
a Ef (z, = Hy) = jz )*p(z) dz ’ 


which measures itS variability about the level of the mean. 
pmmee the probability distribution Oe the same for all 
times t, its form can be discovered by plotting a histogram 
or relative frequency plot of the set of observations 

(Zp rZor-+-rZ,)- Additionally, the mean and variance of 

the stochastic process can be estimated by the sample 


averages taken over time: 


N 
= od 
z= F lL 4 
=] 
and 
N 
2 _ 1 a. 
Ci = > ee) a. 
=i 


Meee BASICS OF DIFFERENCE OPERATORS AND LINEAR FILTER 
MODELS 


This section provides a brief discussion of the differ- 
ence operators currently employed in writing Box-Jenkins 


models. Those most commonly used are described below in 


Table I. 


i 





Operator Name Definition 


BE Backward Shift Operator BZ. ep 
of order m B 
B Z = Zs 
ro Forward Shift Operator f => 
of order m 
ia eer) 
mM = 
: Ze 2+ 4m 
ve Backward Difference da = (l- B)? 
Operator of order D 
and season s ve Tune Bs)P 


Table I. Basic difference operators. 


It 1s noted here that gee could be written as a 


a 
binomial expansion; that is, 
D 
D 7 os AOU, * D ee DES 
Voz, (1 - B) 2, ay) sg) (ee) (—B) )Z, 
yao 
D 
= D) -_,,D-3 
yee 


The stochastic models employed in the Box-Jenkins method 


are based on the idea that a time series in which successive 
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values are highly dependent can be regarded as having been 


generated from a series of independent "shocks" ay: These 


shock terms, sometimes referred to as white noise terms, 

are random drawings from some fixed distribution, generally 
2 ; 

assumed Normal with mean zero and variance oO. The white 


noise process a, 1S supposed transformed to the process z 


t ie 


by what is known as a linear filter, shown in Figure 2 below 


[Ref. 4, p. 8]. The linear filtering operation takes a 
wv (B) 
White Pees Process 
Noise Filter zy 
Ae 


Figure 2. Representation of a time series as 
the output from a linear filter. 


weighted sum of previous observations, such that Zz, can be 


written 


Z = ut a, Us vy Sy + v5 a, reirrc, ica?s 


Woes Ts), an: 


In general, wu is the parameter that determines the "level" 


of the process, and 
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y(B) = 1+,B+ vB : ae 


is the linear operator that transforms a, into z called 


ie ee 
the transfer function of the filter. The sequence 
VyrVor we formed by the weights may be (theoretically) 
finite or infinite; if this sequence is finite, or infinite 
Smemeonverges, the filter is said to be stable and the 


process Zz, Stationary. The parameter uw is then the value 


ie 


about which the process varies. Otherwise, Z 1s called 
nonstationary, and u has no specific meaning except as a 


reference point for the level of the process. 


D. HOMOGENEOUS NON-STATIONARY AND SEASONAL PROCESSES 

In many cases of interest to the time series analyst, 
the process being examined is not Stationary. Instead, 
the probabilistic structure of the process generating the 
time series changes with time. For example, there may be 
a trend or seasonal pattern in the time series. However, if 
the series can be found to exhibit behavior which is somewhat 
homogeneous, then the series can frequently be transformed 
into a stationary series that can be described by Box-Jenkins 
models. 


A series 1S said to be homogeneous nonstationary of 


order dif Wie vz, is a stationary series for some integer 
d. Here V = (1 - B) denotes the difference operator such 
clake ae Va, = Z, 7 2,14- The fact that a given series is 


nonstationary can be determined quickly by examination of 


ae 





the autocorrelation function, to be discussed later. The 
series to be modeled should be differenced until the re- 
sulting series appears stationary or until the procedure 
appears to be making no improvement. 

A time series is called seasonal when it exhibits 
cyclical behavior over time, showing a regular period. 
Examples of seasonal processes include monthly rainfall, 
monthly crop yields, sale volume of bathing suits, livestock 
production rates and energy consumption. Seasonal patterns 
are often easy to spot simply by observing a plot of the 
series, or through knowledge of the process that generated 
the series. However, in many cases where the variability o£ 
the series is large, seasonal patterns will be difficult to 
distinguish from other fluctuations. Recognition of season- 
ality in a Series is important, Since it provides information 
useful in modeling and forecasting; as will be shown, the 
autocorrelation function can assist in recognizing seasonality. 
The Box-Jenkins modeling approach for seasonal nonstationary 
time series of season length s is to first transform the 
series uSing an appropriate seasonal differencing operator 


to “sweep out" the seasonal effect; that is, take 


Meeedete 1S a trend or other types of non-stationarity pre- 


sent, as well as the seasonal effect, the series can be 
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aifferenced again with the ordinary difference operator 
until stationarity is achieved. This transform would be 


written 


The transformed series Ww. is then modeled as a stationary 
series. For seasonal series where complete stationarity 
can not be achieved, a class of seasonal models is avail- 


able; these are discussed later in this Chapter. 


E. GENERAL CLASSES OF BOX-JENKINS MODELS 

This section will provide brief descriptions of several 
of the most commonly used classes of Box-Jenkins models, 
including autoregressive (AR) models, moving average (MA) 
models, mixed and integrated autoregressive-moving average 
(ARMA/ARIMA) models, and a generalized seasonal ARIMA model. 

1. Autoregressive Models 

In autoregressive (AR) models, the current value of 

the series under study is expressed as a linear combination 
of previous series values that explain the current obser- 
vation, plus an unexplained random (white noise) term a,- 
For a stationary series, or one that has been transformed 
to stationarity, the deviation from the mean value for each 


~e 


period, Z = Z, - u, can be modeled as dependent on weighted 


values of the previous deviations from the mean. The AR 


model can be written 





~e 


h = Oa + Paitin + 2.66 + Ans + a, : 

In the expressions which follow the tilde will be dropped, 
but it 1S to be understood thatthe mean has already been 
subtracted from each observation. This equation is called 
an autoregresSive process of order p, where >. represents 
a Scalar weighting coefficient for the hae previous period. 


The model can also be written in more compact form as 


> (B)2, tee Oe 


where 


a polynomial in the operator B. Then, $¢(B) represents a 
general non-seasonal AR operator. The reason these models 
are called autoregressive can be seen from recalling 
regression analysis; here the model relates the dependent 
variable Zi, to a set of explanatory variables (2, 4724297 


me. , Z ), which are previous values of the time Series. 


t-p 
Therefore, the model is autoregressive. For example, when 
the value of the time series at time t depends only on the 
value at time t-l, the model becomes an autoregressive model 


of order 1, denoted AR(1). Similarly, if the value in 


period t depends on the values in periods t-l, t-2 and t-3, 
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the model is AR(3). The following equations are the mathe- 


Matical representation of those two examples: 


st niet 


AR(1).: Zz, = O52,4 eo 


AR (3) +a 


Rice eee os ee a EB t° 


N 


Z- Moving Average Models 


In moving average (MA) models, the assumption is made 
that the current value of the time series at period t can 
be expressed as a linear combination of the previous fore- 
cast errors (or, residuals), a, = Z, - A - The general 


equation for this model can be written 


a eee 


which is called a moving average model of order q, where 
Ig represents the scalar weighting coefficient for the see 
previous period. Like the AR model, the MA model can be 


written in more compact operator notation as follows: 
Zz = 8(B)a, ’ 


1S, 


where 
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a polynomial in the operator B. The MA model implies that 
the analyst can gain valuable information fom future 
predictions by consideration of the weighted sum of a 
number of previous forecast errors or residuals. The 
following are examples of MA models of order 1 and order 


See MA(l) and MA(3): 


MA(1) : Zz, = (l - 9,B)a, oa 954, 4 ; 
MA(3) : = (lsh Boe esos ea 
pet Res igo. 3 t 
Do pa odes 


3. Mixed and Integrated Autoregressive-Moving Average 
Models 


A natural extension of the autoregressive and moving 
average models is to construct a combination of the two. 
Such mixed processes are called autoregressive-moving average 
(ARMA) processes, of orders p and g, often written as 
ARMA(p,q). Box and Jenkins [Ref. 4] notes that for many 
series encountered in practice, the inclusion of both AR 
and MA terms in a model results in fewer total parameters 
than would be required for either a pure AR or pure MA 
process. The ARMA(p,q) model may be written in the usual 


form as 


PAF 





in the more compact operator notation, the same model is 


written 
o(B) Zz, = 9 (B)a, : 
As an example, an ARMA(2,2) model would be written as: 


| 2 
(1- 9,B- $,B*)z, = (1-6,B-6,B)a, , 


or expanded, as: 


+a -=- 8 S) 


Pier iee oc feet oe ti] 2 2t=2 
Recall that with differencing or other transformation 
techniques, a stationary series can be obtained from a 
homogeneous (SeaSonal or non-seasonal) non-stationary series. 
When the differencing technique is integrated directly into 
the ARMA model, the result 1S an integrated autoregressive- 
moving average (ARIMA) model. The model is often written 
ARIMA(p,d,q), where p and q retain their previous meanings, 
and d represents the order of (non-seasonal) differencing. 
The model may also be generalized by including an additional 
constant term oe that will have the effect of adding a 


deterministic polynomial trend of order d. The value of oe 


may be found in general by the formula 
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D> 
Il 

> 

Q) 


where 


ie 
re ) dll = } O51] 
tal j=l 


The general ARIMA(p,d,q) model can be written compactly as 


d = 
(L) o(B) V Zo 85 + 6(B) ay - 
An equivalent form for ARIMA(p,d,q) models that is 


sometimes seen is called undifferenced form: 
¥(B) z eg stot SE 
where 
ns) ee ee) (Es =e (0) 412) Oe. s)) Sa 


The "differenced" form, equation (1), is the form usually 
seen. In the differenced form, the transfer function 9(B) 
is assumed to be stable; that is, all the roots of the 
difference equation $(B) = 0 are outside the unit circle. 
Clearly, ¥(B) is not a transfer function of a stationary 
series when d > 0, since ¥(B) has dad roots on the unit 


errcle. 
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4. Seasonal ARIMA Models 

The general class of ARIMA models previously dis- 
cussed can be further generalized to allow for the modeling 
of seasonal models with the addition of appropriate seasonal 
operators and parameters. The seasonal model is called 
an ARIMA model of order (p,d,q) x (P,D,Q) ., where p,d, and 
Geretain their previous meanings, and P, D, Q and Ss are as 
defined in Section A of this chapter, page 17. The ARIMA 
(p,d,q) x (P,D,Q). model can be written as the product of 


nonseasonal and seasonal operators: 


d 


Sie) = S 
o(B)(B ye Z, 0 = oA spans) (Us) 1S) (0s) ja, 


As an example, consider the ARIMA (1,2,1) x (1,2,1),, model: 


(1-$,B) (1-8 elié\actys 2 seg. & (1-8 ,B) (1-0 Bl)a 


ile 1s O ieeZ t 


IE THE AUTOCORRELATION AND PARTIAL AUTOCORRELATION FUNCTIONS 
Mathematically, the autocorrelation at lag k, denoted 


Pre is defined as: 


The autocorrelation function of a time series describes the 


association (mutual dependence) among values of the same 
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series taken at different time periods, the difference in 
time being referred to as the lag, denoted k. Autocorrela- 
tion of itself implies nothing about a change in one 
variable cauSing a change in another. However, the auto- 
correlations may provide important information about the 
structure of a data set and its pattern. Ina set of com- 
pletely random data the theoretical autocorrelations among 
Successive values will always be zero, whereas data values 
of strong seasonal or cyclical behavior will be highly 
autocorrelated. MThis P, can be estimated using the time- 


averaged sample autocorrelation, denoted r defined as: 


eu 


N-k _ _ 
) (2, - 2) (z - z) 
t=1 


L 
N t+k 


ba ae ~ 
k 57 
Z 


A plot of the autocorrelation function versus the lag k, 
called a correlogram, is very useful for the purpose of 
determining if a process is stationary and for identifying 
the appropriate model. 

Another function useful in identifying the appropriate 
model for a given series is the partial autocorrelation 
function, often called the pauto function. Partial auto- 
correlations are analogous to autocorrelations in that they 
indicate the relationship of the values of a time series to 
Various time-lagged values of the same series. However, 


they differ from autocorrelations in that they measure the 


ee! 





strength of the relationship between values of the series 
of various lags after the effects of other lags have been 
removed. In effect, they show the relative strength of the 
relationship that exists for varying time lags. For time 
series models of the types to be considered, the partial 
autocorrelation coefficients can be calculated several 
ways. One method of calculation, described by Box and 
Jenkins [Ref. 4, Sec. 3.2.5] begins rae the following 
equation, satisfied by the autocorrelation function, where 


the OL are autoregressive parameters: 


ie kieg-2 Ok Haj tr © Pek Pg-e 6 IT Freee: 


This leads to what are known as the Yule-Walker equations, 


which may be written 


> cco oan Pres Py 
ep, 2 ee > Pr9 Po 
0 0 0 1 t 9 

eon pee oe es ee k 


Seuvaing these equations for k = 1,2,3, ..-, successively, 


we obtain 


SZ 





b “aL 
2 
a 2 SO 7a 
22 Ze ie 
ik Py 1 - py 
Py IL 
1 Py Py 
4 iE 05 
0 0 0 
2 1 3 
o = ee orc. 
32 ty 0. 
, 2 i 
25 5) 1 


In general, for Pret the determinant in the numerator has 

the same elements as that in the denominator, but with the 
last column replaced with Py The quantity PK? a funcelon 
@emcne lag k, is called the partial autocorrelation function. 
Other methods for calculating partial autocorrelations include 
Successively fitting autoregressive processes of orders 

1, 2, ...-, k by least squares and picking off estimates of 
Dose De Bek - by of the last coefficient fitted at each 
stage. Another method, discussed in Appendix A3.2 of 

Ref. 4 and due to Durbin, generates estimates of an auto- 

' regressive meocese of order k+l recursively from estimates 
of autoregressive processed orders k and less. It is 

derived by observing recursive relationships from the 


Yule-Walker equations. The recursive formulas are: 
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nw 


Ol, 4 k. 2 Piel, kt1 9k, kt : 


fie, 2, ..., K, and 


k 
Rn "k+1 - es *k+1-3 
Pe+1, k+l ; * 
i- Ome Io 
j=1 hy 


Another way of looking at partial autocorrelations 


mento consider a time Series {z,} and the model 


Z. <5 la) ee Leeer/ + b.Z a ees Hees ees hee ; 


where: the b values are the least Squares estimates of the 


linear regression coefficients (8's) in the model 


Zz = Se os 


: : By 2,4 + ... ¢ 


oo SGee Ty ae 


Let ae be the residual of z, after removing the linear 


= 


effect of z fOr. Z Such that zZ = Z - 2Z 


cena 2) eles t! t & t 
Then, the partial autocorrelation of lag k, denoted Tae ee 
is defined to be the simple autocorrelation of lag k for 


~e 


the adjusted series (Zi 1Zor Sane: Z,)° 
The exact expression for partial autocorrelations of 
a moving average process is complicated, but an approximate 


duality between autocorrelations of an autoregressive 
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process and pautos of a moving average process implies 
that the pautos of a moving average process would have the 
same general behavior as the autocorrelations of an auto- 
regressive process of the same order. This approximate 
duality also implies that the pautos of an autoregressive 
process of order p will have the same general behavior as 
the autocorrelations of a moving average process of order 
p.- Thus, one function can be examined and used to confirm 
the other, and may be of use in model identification, discussed 
in section G of this chapter. 
fee COMMENTS ON MODEL SELECTION EMPLOYING AUTOCORRELATIONS 

AND PARTIAL AUTOCORRELATIONS 

The purpose of this section is to provide a brief over- 
view of the methodology for selecting the proper model for 
a given time series, through both examination or prior 
knowledge of the data itself and examination of plots of 
the autocorrelation and partial autocorrelation functions. 

A not yet stationary series should first be differenced 
or transformed until it appears to be stationary; character- 
istic of the autocorrelation plot (correlogram) of a non- 
Stationary series is a very slow damping out of the auto- 
correlations. When this property of the correlogram is 
observed, the analyst should difference the series until no 
further improvement is reached. For ARMA(p,q) models, the 
tentative identification of the model class, that is, the 


determination of p and q, can be done by examining the 
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sample autocorrelation and sample partial autocorrelation 
Mij@ertons OL the given time series with the theoretical auto 
and pauto functions of members of the general linear class. 
For most Stationary time series, an adequate fit can be 
found in a model with p and q relatively small, say three 


or less. 

The following figure displays examples of autocorrela- 
tion and partial autocorrelation Tegue aS plots for several 
basic model types, described in Ref. 17. 


Figure 3. Examples of Auto and Pauto Plots 


Figure 3(a). Auto and Pauto vs. Lag for AR(1) Model 
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migurers(D)-  AutG and Pauto vs. Lag for AR(2) Model 
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Figure 3(c). Auto and Pauto vs. Lag for MA(1) Model 
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Figure 3(d). Auto and Pauto vs. Lag for MA(2) Model 
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Figure 3(e). Example Auto and Pauto vs. Lag (ens 
Mixed ARMA(1,1) Model 
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As can be seen from the figures, there is a wide variety 
of possibilities for appearance of the correlogram and 
partial autocorrelation plots for even the simplest of 
models. Sometimes the pattern of estimated autos and pautos 
can be quickly classified in terms of one of the simpler 
basic models; however, most real data will generate auto 
and pauto plots that resemble those of Figure 3(e), those 
of a mixed model. Here some free association may be re- 
quired to infer a pattern from the correlogram, or more 
than one pattern may be implied. However, most of the time 
a selection of something 1s possible. As is evident, this 
is a somewhat subjective process, where the quality of 
identification will improve with experience. 

The following table, adapted from Box and Jenkins 
[Ref. 4], summarizes the properties of autoregressive, 
moving average and mixed ARMA processes; an understanding 
of its contents will be helpful as thumbrules for initial 
model identification. 

It can be noted from examination of Table II that there 
exists what can be termed a "duality" relationship between 
the autocorrelations and partial autocorrelations of pure 
autoregressive and moving average processes. For example, 
the plot of autocorrelation versus lag for an AR(p) process 
would appear the same as the plot of partial autocorrela- 
tion versus lag for an MA(q) process, for p = q. Similarly, 


the plot of partial autocorrelation versus lag for the AR(p) 
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Autoregressive 


Process 

Model in terms _ 

of previous o(B)Z, = a, 

Zz, values 

a 

Autocorrelation Infinite and tails 

Function off; composed of 

damped exponentials 
and/or damped sine 
waves 

Partial finite; that is, 

Autocorrelation there will be p 

Function non-zero partial 
autocorrelations 

maole IL. 


Moving Average 


Process 


‘ia = 8 (B)a, 


finite; that 

is, there will 
be q non-zero 
autocorrelations 


Infinite and 
coats OLe- 
dominated by 
damped exponen- 
tials and/or 
Sine waves 


Mixed 
Process 


o(B)Z, = 6(B)z, 


Infinite and 
tails off; 
composed of 
damped exponen- 
tials and/or 
damped sine waves 
after the first 
q-p lags 


Infinite and 
tails off; 
composed of 
damped exponen- 
tials and/or 
Sine waves after 
the first p-q 
lags 


Properties of Basic ARMA Processes 


process would appear the same as the plot of autocorrelation 


versus lag for an MA(q) 


process, where p = q. 


For seasonal series, where differencing has been unable 


to remove all the nonstationary seasonal behavior, 


the auto- 


correlation and partial autocorrelation plots will generally 


exhibit a "spike" at lags equal to integer multiples of the 


period of seasonality. 


For example, 


with a period of 12, then the autocorrelations 
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1f the data are seasonal 





co Lag cee and the partial autocorrelations 


nD” Doe a 


o42,12° oo4,24' 036,36" »«» would be amplified. This would 


indicate a need for one or more seasonal parameters to be 
included in the model. The following figure is an illus- 


tration of such an auto/pauto plot. 


Figure 4. Example Auto and Pauto versus lag for mixed 
ARMA (p,q), with seasonality remaining of 
period s = 6. 
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H. ESTIMATION OF MODEL PARAMETERS 
Suppose a time series has been tentatively identified 


as an ARIMA(p,d,q) model: 
@(B) Vz = §@ + 6(B) a 


There are then p+q+2 unknown parameters of the model to be 


meammated, (>,.,0,, .-.+-, 6.,9_,9.,8 eee waits). The 
lea ba Oo 1 a 


Phd q 
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Box-Jenkins procedure separates this estimation into two 
parts. First, estimates are obtained for the AR and MA 
parameters d= (bd, 0b5: - sony >) and 6 = (85755, emo.) > 
then estimates are calculated for 8 and of, which are func- 
jeons of 6 and 6. 

The uSual procedure is to select those parameter values 
6 and § that minimize the sum of squared model errors 
(residuals). Let ea 8/0 ~ os - b> ee es. = ,,) and 


~ 


_ ae = ; _ 
Ww, = V (Z, Wo): Then, it can be shown that Bee EW 3 


and that the model can be rewritten as 


$(B) (Wo - uw) = 8(B) a. 


OL 


-1 . 
(B) o(B) (W, - u) 


a 
I 
D 


Now, let W, the sample mean, be the estimate of Bae 


mee) Oo, then We 1s usually zero). Also, let 
a, = 6 +(B) $(B) (W, - W). 
ie tc 
Set 
N 
a(6,0) = yen 
ll t ! 
t=1 


where N = n-d is the length of the differenced series. The 
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objective is to select those parameters ¢$ and 8 such that 
S($,8) is minimized over all values of the (p+q) -dimensioned 
parameter space. Since the equation in S is nonlinear 

in the parameters, iterative search methods are generally 
used to determine optimal model parameters. Having found 
these optimal parameters and the minimum S, of can be 


estimated using 
of = $§(6,6)/N-p-q . 


Chapter III contains a more detailed description of a 
nonlinear algorithm to estimate parameters of a generalized 


seasonal or nonseasonal ARIMA model. 


meee LAGNOSTIC CHECKING OF MODELS 

After the model has been tentatively identified and 
parameter estimates have been calculated, the next task is 
to test whether or not the original model specification 
was correct and that the model itself is adequate in fore- 
casting power. The process of testing the model can take 
many forms, but will usually include at least the following 
two steps: 

1. Generate a simulated series from the estimated 
model and compare the simulated series and its 
auto and pauto functions with the original series 
and its respective auto and pauto functions. This 


comparison is essentially subjective. 
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2. Calculate the residuals of the estimated model, the 

a.'S, and compare the properties of the residuals 
with those assumed for the shock terms of the actual 
process. The residuals should be normally ee 
buted and uncorrelated with each other; that is, there 
should be no discernable structure in the residuals. 
There are many quantitative statistical tests and 
data analytic tools that can Be eo1ie to the 
residuals to test hypotheses of normality and zero 
autocorrelation. 

A plot of the autocorrelation and partial autocorrelation 
functions of the residuals will provide not only a test of 
whether or not the residuals are uncorrelated, but, if they 
meee cOrrelated, the plots can be used to suggest improve- 
ments to the model. For example, suppose the model was 


tentatively specified as the ARMA(1,1) model below: 


Suppose also that the autos and pautos of the model residuals 


. = UL, where the u, S 


os : 2 
are white noise, uncorrelated with variance ae Then, 


Suggested the model (1 - 0.3B) a 


these two models can be combined into the ARMA(2,1) model: 


Salo eter On55) (20 — 2)) = {1 + 0.7B) u 


1e t 


This should be a refinement over the original model. 
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J. FORECASTING 

The objective in forecasting is to predict future 
values of the time Series with as little error as possible. 
The criterion most often used for forecasting 1s to calcu- 
late that forecast which ailaiaieec the expected mean square 
forecast error. Therefore, if 2, (2) denotes the forecast 
for lead time 2 from time origin t of the value Zags 
the objective is to find zie) such that the objective 
fonet 1On 


A Z 
EC(Z, 45 ~ Zz, (2) ] 


1S a minimum. This forecast is given by taking Z (2) as 


Z 


given Zyr seen 2h! 


the conditional expectation of Zi yge 


N > 


CS =e OE 


Zeng! Z E72 E a1! 1 
The forecast can be easily generated recursively from the 


mathematical model, utilizing the facts that 


EZ, _s|Zp0-+++2y] = ee : 
momejy = 0,1,2... (where t is the current time), and 
EZ ys lagers Zz] Sen cor ge= 1,2; 
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and 
0 Eom tlmes a t 
a ee ENE Se 
For example, consider the Boe euets model 


iL oo Ole Wee z. = (1 + 0.3B) a, 


which could be written in the expanded form as 


Z. = Q.52 ~ 0.6z,_. ay + 0.3 ayy x 


Now, assume it is known that z Ss ee sel .0, and 


100 “96 


the calculation has been made to obtain the residual 


4100 = 0.2. Then the forecasts of Z5o1! 2492 and 2493 


made from an origin of t = 100 can be calculated from the 
model as follows: 
mics, tor § = 1, set 


nN 


Zigq(t) = ElZy 9121 99"%99"°7 7741!» 


EL0.52, 49 - 0.624, ~ As01 ~ 0.38549! j 


% 


N> 
- 
iI 


100 


and 
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top (2) == (0-5) (1-4) = (0.6) (1.0) + 0 + (0.3) (0.2) ; 


then our forecast 
Zan) =—— 057 =eOno ImOs06m = O2160.; 


Next, for 2 = 2, set 


aw 


Z199(2) = ElZj991%191"21007299" °° 7121! » 


El0.52Z, 5) - 0-629 + A102 7 0.385) 


21909 4) 


ore . 
en(2) = (Ge) Waa ayy CD = (005) Carag a Ol ae Oe 
Ai ae) = (0.5) (0.16) - (0.6)(1.4) = -0.76 ; 
meainitlarly, for & = 3, set 
21993) = El219312192'2101'7100,2997 "7/21! 
ara) = EID 5Z 5 45 - ORGz yaa Sere Oe ania 
121 99'2%997°° 7124! 
Bans) 9 = (GES (Za ia 2) ) OG) (25 ayy CED aCe 
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and then 


An 


Z4 90 (3) ccm) 0. 7/6) )—- (020) (0.16) = Ole ) 


Clearly, this process can be continued into the future as 
long as deSired, recognizing that the expected forecast 
accuracy will decrease as 2 increases. 


Let 
e, (2) = (2 = 2 Guy) 


denote the forecast error for lead time 2% beyond the forecast 
time origin t. It can be shown that e, (2) is given by 


the relation 


e, (2) = a + 


pea) WERE SG) ee ee 


where the scalar weights Ds are determined from the equation 
IL =x 


Tee Oe (BieGLo— B) — 6 (B) 


uSing the known parameter values. The variance of the 


forecast error is given by 


2 _ Z 2 a ae ae 
meee) ) = (1 +p, +... +H) 4) (05) = als l URion 
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From these relationships, a confidence interval of n 
standard deviations about a forecast of lead time % would 


be given by: 


ve | 1/2 


ci (2, (2)) = 2,(2) + nfl + J v1 ee 


jal 


The value of n could also be though of as the percentage 
point for the desired level of confidence using the Normal 
distribution [Ref. 4]. It can be noted from the definition 
of e, (2) above that the one-step-ahead forecast error, 


e, (1) is simply the value a that is, 


t+’ 


This is an explanation for the use of the term "residual" 
to refer to the white noise or shock terms. Additionally, 
from the foregoing it is evident that the forecast error 
Variance is a non-decreasing function of the length of 
the forecast lead time %; therefore, the confidence bands 


must become wider as the forecast lead time increases. 
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III. GENERALIZED ARIMA MODEL PARAMETER ESTIMATION 


The purpose of this chapter is to describe the compu- 
tational details of the algorithm WMARQRDT, employed with 
appropriate subroutines to estimate Box-Jenkins parameters 
of a generalized seasonal (or non-seasonal) ARIMA model. 


This generalized model can be written parsimoniously as 


follows: 
Ss _ s 
6(B)$(B )W, = 8 +8 (B) O(B ya, ; 
where 
oa Dee 
mee is V Vqzt! 
and 


Zp is the transformed time series (for example, 


_—— oo 
Zs ln Z, Or Z, v2.) ; 
The model employs a Marquardt-type [see Ref. 13] non- 

linear least squares algorithm for determination of the 
model parameters, searching in the parameter space for the 
set of parameters which minimizes the sum of the squared model 


residuals. 
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A. MATHEMATICAL DISCUSSION OF THE MARQUARDT ALGORITHM 
Reference 13 describes Marquardt's methodology for 
development of a general algorithm for Pee er scares esti- 
mation of nonlinear parameters. Most algorithms for least- 

Squares estimation of nonlinear parameters have been cen- 
tered about either a pure linear iterative model based on 

a Taylor series expansion or using some form of the method 
of steepest descent (maximum negative gradient). Since 

both methods have severe potential pitfalls (Taylor series 
due to divergence of the iterates due to an unfortunate 
choice of initial values, the steepest descent due to slow 
convergence after the first few iterations), a method termed 
"Maximum neighborhood" method was developed by Marquardt 
which, in effect, performs an optimum interpolation between 
the Taylor series method and the method of steepest descent, 
the interpolation being based upon the maximum neighborhood 
in which the truncated Taylor series gives an adequate 
representation of the nonlinear model. 

As discussed in Ref. 14, the idea of Marquardt's method 
can be briefly explained as follows. Assume the algorithm 
begins with a vector of initial parameters ero aren) 
If the method of steepest descent is applied, a certain 
vector direction Rae where the subscript g represents the 
gradient (which gives the direction in which the rate of 
change is the greatest), is obtained for movement away from 
the initial point values in the parameter space. Due to 


potential nonlinearities in the sum-of-squared-residuals 
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surface, S(6579019679,) 5 in the parameter space, Re may 
Bemene best local direction in which to move toward opti- 
Pedity (minimal sum of squared model residuals), but may 
not be the best direction. However, the best direction 
must be within 90 degrees [Ref. 8] of fo or else the 
value of S(o,%,98,9) will increase locally. Now, the 
linearization, or Taylor series method, may result ina 


different correction vector, h Marquardt found experi- 


b° 
mentally that for a number of practical problems, the angle, 
say €, between the vectors a and Ay edi ieathc piscine 

80 degrees < § < 90 degrees. In other words, the two 
directions were nearly normal in most cases. The Marquardt 
algorithm then, provides a method for interpolation between 
the vectors as and he, employing a sutiable step size, which, 
in general, converges more quickly than either of the two 
“parent" methods alone. The remainder of this section des- 


cribes the mechanics of the algorithm, as applied to the 


estimation of Box-Jenkins parameters. 


B. MODEL INPUT INFORMATION 
The following is the required input information to 


the WMARORDT parameter estimation model: 


Dw et Liew ratnstormed time series (Lor example, 
i] = | = 
WwW, in Zi, OL We v2, 
N the length of the time series 
S the length of the seasonal period (for no 
seasonality, s = 1) 


S/S) 





P the number of seasonal AR parameters, P > 0 


p the number of non-seasonal AR parameters, 
p> 0 
Q the number of seasonal MA parameters, Q > 0 
g the number of non-seasonal MA parameters, 
q > 0 
NDIFNS the number of non-seasonal differences to 
be taken 
NDIFS the number of seasonal differences to be taken 
$ = (® pee © ) initial estimates of seasonal AR 
Ss P,s 
parameters 
$ (by rbogreee rd) initial estimates cf non-seasonal 


AR parameters 


6 = (0, Eo. avi: 96 = initial estimates of seasonal 
f f é MA parameters 
6 = (6,,8,, cae 6 ) initial estimate of non-seasonal 
q MA parameters 


C. ALGORITHM DESCRIPTION 

Initially, the model accepts inputs of the time series 
length and the (transformed, but undifferenced) time series 
itself from the user's disc, where it is stored in File 
FTO2FOQO1. The remainder of the input data is then entered 
interactively on the terminal as prompted by the model. The 
model then forms the initial parameter estimates into vectors 
of lengths p, P, q, and Q, and employs them to calculate 


an initial model residual sum of squares. 
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Prior to initiating the sum of squares calculation pro- 
cess, the model unravels the polynomial of AR and MA operators, 
and forms a vector of coefficients for both the AR and MA 
Sides of the model. For example, consider the relatively 


simple (2,1) x (2,1),, ARIMA model, written 


2 1 epee IZ 
(1-9, B-$5B°) (1 - , 49B°)W, = (1-8,B-8,B°) (1-0) 1B “ay. 
Unraveling, the model becomes first: 
2 _.. eee 
Ree oe M8) yop) = (2-8, B-8,B ) ta.) 12%e-12! 5 


and then takes the form: 


W,_-® 


Mm 12 ee 1, eH 13 econ oop aa esia 


9 850 


I ee ee ier 


Bese cee 1 eal 112 t=13° 
This 1s represented in the computer as four vectors: two 
vectors, mt and y, composed of the expansion AR and MA 
coefficients respectively, and two vectors composed of the 
indices of the time series terms associated with the coeffi- 
cients. For the model in the example above, these vectors 


would be written as follows: 


INDEX... =eOrmle2 ai 1359 2, 14) 
Y = (118) 4970, 7 (85) (Oy 49) 1950 (85HG 45)) 


S) 5) 





INDEX, —mertowme tl, iit3, 2, 14). 

This format for unraveling of a generalized product of 
polynomials of seasonal and non-seasonal parameters permits 
Simplified calculations in the determination of the sum of 
Squared residuals, since the arithmetic need only be per- 
formed for the non-zero parameter values, those contained 
in the T and y vectors. In the actual calculations, per- 
formed in FORTRAN, the index vector origin is translated 
to the right by the amount p + Ps + q + Qs, called IADDIT 
in the model, in order to prevent obtaining zero or negative 
subscripts. | 

The calculation of the residual sum of squares is a 
three stage process. First, the model calculates ten 
values of e,, the "forward" white noise terms, uSing the 


te 


following relationship: 


The assumption is made that (e grees rey) are all equal 


lolpaeal 


to zero. Having found the desired values for (2,1 egr ee ey) 


we can now backcast the necessary values of W in order 


t’ 
to enable us to calculate the required white noise terms 


a as follows: 


t’ 
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pt+Ps gq+Qs 


ea ) heer 
i=2 j=2 


solving for (Wo, W_y/W Now, with the 


-2°°°° 4) _ (p4sp+q4sQ))° 


W's backcast, we can proceed to calculate the estimated 


white noise terms (residuals) uSing this relationship: 


pt+Ps qtOs 
ese TW t) ) Y 594-3 
i=2 4=2 


The values of ay are calculated for t = l, ..., N. Once the 


a,'s are calculated, it is a simple matter to calculate the 


sum of Squares, uSing 


2 


N 
S(¢,%,69,0) = ) a. “ 


t=1 


Having calculated the reSidual sum of squares for the 
initial parameter estimates, we form the parameters into 


a vector 8, such that 


where k = p+P+q+Q ; that is, 
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The subroutine PARSH then calculates the derivatives 


Ja 
“i,t ~ OB ’ 
i 
over all values of t = 1,...,N and 1=1,...,k. Using the 


residuals aren calculated earlier, the derivatives are 
estimated numerically uSing a perturbation of each parameter 


of amount 6, such that 


“it = ta, (Byr---7B,)-a, (By 7--+ 78; t5,--- 18) T/6 


Then with ta, f and x. Supplied for the current parameter 


rt 


values, the following quantities are formed: 


1. The kxk matrix A = {aj5} , where 
N 
a ) (Xi 4) Oe) 
t=1 
Meeeethe vector G = tg, } , where 
N 
a, = Tay) fay 
t=1 
3. the scaling quantities D. = vA. : 
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Then, the modified (scaled and constrained) linearized 


equations 


aia) = Sig] 
are formed, where 
si” = 5 af De , LF ; ; 
A, ,* = A, ,* + € &= .01 , 
and 
Sree Shae 


It can be noted here that the Gog a Mattei NS) eine Cite wel 
correlation matrix for the model parameters. Those equa- 
tions are then solved for h*, which is scaled back to give 


the parameter correction vector has where 


h. = h.* Ds, e 
J J / y 


Then, the new parameter values are formed from the h vector, 


men that 


Se, 





and the sum of squared residuals for the doe Ss he 
parameters is calculated. 


Now, if S (Byarg) < S(B ), the parameter corrections h 


OLD 
are tested. If all elements of the h vector are smaller 
than some epSilon, say .0001, then convergence is assumed 
and the kxk matrix fax] is used to calculate the covari- 
ance matrix of the estimates; otherwise, Bop is reset to 


BNEW? the constant € is reduced by a factor F, say .l, and 
computation returns to the calculation of a new set of 
derivatives. 


However, 1f£ S(8 ), the constraint parameter 


new) 7 3 8orp 
—€ 1s increased by the factor F, and computation is 
returned to the stage where the A matrix is formed. In 
all but unusual cases, a minimal sum of squares will be 
found. However, an upper bound can be placed on €, and 
when it is exceeded, the search 1s terminated. When con- 


vergence is reached, the desired output information is 


calculated. 


D. MODEL OUTPUT AND DIAGNOSTIC CHECKS 


The output includes: 


me B= (Byre-+1 By) the least-squares estimates of 


Qyrere rd 
® 


D 


1p 


po 
q 


mo 


Gi reee 


C Pac 


1908 0s 
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2. co. the residual variance, Toe a0) S(¢,9,98,9) 
BS. V the covariance matrix of the parameter 
estimates, formed from 
= 7 WU Sih SV 
4. S; the standard errors of the parameter 
estimates, 
=e = aie r a = Je Be are creel Ors 
D Ri the correlation matrix, obtained from 
Re ey ee. 
1) 13 il Jj) 
6. g the overall constant term, where 


O ia a 
oS =u G =, for 


p P 
G = (1- J 6,)(1- J} $5) ana 
1=1 j=l 
N 
ueaye| 
ae aoe 
=1 


Finally, diagnostic checks on the model are performed. 
This includes calculation of residual autocorrelations 


mea\k), obtained from these formulas: 


1 lS) ae SOS) otars 1 09)) - 


Ger 





where 


and k goes from 1 to the maximum desired lag K, usually 
40. Additionally, a chi-square statistic is calculated 


from 


K 
= Zn 
=f ee as) ; 
k=1 


2 
Xv) 


and 1S compared with a chi-square distribution with 
meee DP - G- P- Q degrees of freedom. 

Examples of input methodology and appearance of output 
for this model are contained in later chapters of this 


paper. 
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nye PBSC RIETION OF THE TIME SERIES EDITOR 


In this chapter, descriptions are given for each of 
those programs in the Time Series Editor which interact 
with the user. There are fifteen separate program modules 
currently included in the Editor; there are other programs 
utilized by the Editor that are completely transparent 
to the user, and are not described here. However, listings 
for all programs contained in the Editor are provided in 
Appendix D to this report. 

No attempt 1S made in this chapter to describe the 
actual mathematical calculations or algorithms that are 
performed by the programs. Rather, the objective is to 
provide general descriptions of what each program can do 
for the user, and how the user interacts with the programs. 
The programs described in this chapter are, in the order 
Zeesented: TIMESER EXEC, ZFORMAT, CMSWORK, TRANS, DIFF, 
PLOT, AUTO, ESTIMATE, YESTSEAS, WMARQRDT, XSUMSQ, FORECAST, 


ROOTS, SIMULATE and GENERATE, and HELP. 


A. THE EXECUTIVE PROGRAM 

The heart of the Time Series Editor is a master program 
Seebed TIMESER EXEC which provides file control for all of 
the other program modules, controls input and output, super- 
Vises the necessary CP/CMS protocol and provides instructions 


interactively to the user concerning the contents of the 


63 





Time Series Editor and how each program can be used. 
The TIMESER EXEC program is written in the IBM-360 CP/CMS 
Executive language. It 1s the only program in the package 
not written in FORTRAN; consequently, it should be the only 
program that would need modification if the Editor were to 
be adapted to another FORTRAN-capable time-sharing system. 
After the user has logged into CP/CMS and linked to 
the disc space containing the Time Series Editor (instruc- 
tions for this procedure are provided in detail in the 
User's Guide to the Time Series Editor, included as Appendix 
A to this report), the entire Editor package is made avail- 
able to the user by the command TIMESER. On entry of this 
comand, the executive routine, TIMESER EXEC, can provide 
a guided tour through the Editor. It will briefly des- 
cribe what the Editor can provide and asks the user what 
Option he wants to use. On the basis of the user's response, 
the Editor then advises the user what input data is required 
and how it should be entered, either through a data file 
entered offline or parametric data entered interactively 
at the terminal. When the user selectS an option for exe- 
cution, the TIMESER EXEC routine loads the appropriate pro- 
gram package and automatically manipulates the required 
input and output data files. For the more experienced user 
who does not require detailed user instructions, there 1s 
a shortened version of the TIMESER EXEC. Upon entry into 


the shortened version, the user iS immediately asked which 
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option he desires. This version of the Editor is entered 
by adding the argument "s" when logging into the Editor; 
that is, type TIMESER S when logging in. The TIMESER S 
(shortened) version of the Editor provides exactly the 
same program options to the user as the longer TIMESER 


version. 


B. DATA AND FILE MANAGEMENT 

Whenever data is required, the user is prompted at 
the terminal by either the TIMESER EXEC or the program 
module being executed. In most cases, the necessary user 
response is a short alphanumeric character input during 
execution uSing the terminal keyboard. However, when the 
time series data itself is entered, the user may enter the 
data by card deck offline. Similarly, most output is pro- 
vided right at the user's terminal. Output such as listings 
of transformed series or plots are sent to the offline 
printer to conserve time and provide hard copies of the 
results. Detailed descriptions of the input and output 
requirements of each program module are provided in this 
chapter, and also in abbreviated tabular form in Appendix 
A, the User's Guide. However, there are some general prin- 
ciples of data management that apply to all programs. 
These are described in this section. 

fee OLttline Data Entry 

The Time Series Editor requires that the user's 


time reside in FILE FTO2F00l1. The data can be read offline 
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via cards and directed to FILE FTO2F001. An example of the 
required format forthe input card deck for the Time Series 
Editor is provided in Figure 5. It is noted here that the 
data format expected by all programs is FORMAT(5F15.6), 
Meeeeded by 4a card containing the length of the time series 
in FORMAT(I3). 
2. ZFORMAT Program 

In many cases, the user will want to analyze time 
series data that are already available ina format other 
than that required by the Editor. To spare the user the 
tedious task of retyping his deck in the required format, 
the ZFORMAT program was written. ee program converts a 
time series data file (FILE FT03001) written in any FORTRAN 
format into the proper FORMAT(5F15.6) onto FILE FTO2FO001. 
The length of the series is entered aS an input on the 
terminal. The data are then ready to use in the Editor, 
With the original file left in FILE FTO3F001. Figure 6 
Shows a sample input deck for such non-standard data. 


Figure 7 gives a sample session where ZFORMAT is used. 
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timeser s 


ByveR LETTER FOR OPTION YOU WANT. 
Z 
IS YOUR DATA IN FILE FTO3F001? 


a 
feoCUTION BEGINS... 


moe YOUR TIME SERIES DATA NOW IN FILE FTO3FO001? 
y 


fveek THE LENGTH OF YOUR TIME SERIES VIA FORMAT 13. 
010 


NOW ENTER THE FORTRAN FORMAT FOR YOUR TIME SERIES DATA, 
PmIemhUDING PARENTHESES; FOR EXAMPLE, TYPE: (5F15.7). 
meee, 2 


weokR FORMAT IS: 
FORMATF10.2 

Boe luls CORRECT? 
n 


NOW ENTER THE FORTRAN FORMAT FOR YOUR TIME SERIES DATA, 
miesJOING PARENTHESES; FOR EXAMPLE, TYPE: (5F15.7). 
eed. 2) 


moUR FORMAT IS: 
PORMAT (FLO. 2} 

Moet diis CORRECT? 
of 


meok DATA IS NOW TRANSFORMED INTO THE PROPER FORMAT 
eer USE IN THe TIME SERIES EDITOR,LOCATED IN FILE 
mmeezFOOl. YOUR ORIGINAL DATA IS IN FILE FTO3FOOIL. 


DO YOU wANT TO GO AGAIN? 
n 


CONTROL RETURNED TO CMS 
R; 


Figure 7. Sample user session with ZFORMAT program 
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3. Interactive Keyboard Input 


a. Numeric Input 
When the user is prompted at the terminal to 

enter such numerical values as the number of AR parameters, 
estimates of parameters or the number of differences to be 
taken, he must enter these values before program execution 
can continue. The user should enter the data according to 
standard FORTRAN practice. For example, integer data should 
be entered (without a decimal point) for counts, such as 
the length of a seasonal period, and for names beginning 
with the letters I through N; floating-point data, such as 
initial estimates of parameters, should be entered with a 
decimal point. Because . typed decimal point overrides a 
floating-point format, it 1S not necessary for the user 
to concern himself with the format for floating-point input 
data. However, care must be taken when entering integer 
data because it must be right-justified in its format field. 
The user is advised in each case when the format is other 
than Il. For example, suppose the user desires to perform 
analysis on a series having only 90 observations. If the 
program that he is executing requires the length of the time 


series, the user will receive the following request: 
ENTER LENGTH OF THE TIME SERIES VIA FORMAT I3 


The user would then enter: 


rae 





column 23 


b9Q (where b represents a blank space). 


Tf the blank were omitted, the program would read the 
series length as 900, and problems in execution would 
Sccur. 
b. Alphabetic Input 

The Editor often asks the user to respond with 
alphabetic input. Generally, this is in the form of the 
response to a question requiring a yeS or no answer, or a 
ie pense to the question of which program to execute. The 
editor has been programmed to read only the first letter 
of such responses. Therefore, the user need only enter the 
first letter of the response for each such inquiry. For 
example, the user should enter Y for yes, N for no, P for 
the PLOT option, X for the XSUMSQ option, W for the 
WMARQRDT option, etc. The other alphabetic input is 
generally entry of titles to plots; these can be any 
combination of numeric and alphabetic characters, as long 
as they do not exceed 72 columns in length. For example, a 


Suitable title would be: 


PLOT OF AUTOS AND PAUTOS FOR SERIES C DATA (2 ORDINARY DIFFERENCES). 


ah Output 


Most of the results are written out at the user's 


terminal. However, in some cases the output is printed 
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offline to conserve time. Such results as plots and 
transformed time series are written onto various files 
(FTO3FO001, FTOSFOO1L, FTO8FO01), and the TIMESER EXEC pro- 
gram causes them to be automatically printed offline under 
the user's identification number (USERID). The plots and 
files can also be printed out at the terminal using the 
usual CMS commands in the CMS mode, outside the Editor, or 
in the CMSWORK program (see next Bae conve However, the 
user should be aware that plotting graphs at the terminal 
is generally ane a slow process. 
5. CMSWORK Program 

This program permits the user to enter the CMS 
environment to perform routine CMS functions without leaving 
the Time Series Editor. Commonly used CMS commands include 
File name alteration, file erasure, file offline printing 
or punching, obtaining a disc status, and listing all disc- 
resident files. However, essentially all CMS commands can 
be issued while in CMSWORK. Figure 8 provides an example 


of the use of the CMSWORK program. 


C. TRANS PROGRAM 

The TRANS program takes a given time series in File 
FIO2FO001 and performs a transformation upon each data ele- 
ment. The options include a change of scale, a power trans- 
formation, a shift of the origin, and a natural logarithm 


transformation. The original data is written onto a file 
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timeser s 


ENTER LETTER FOR OPTION YOU WANT. 
c 
ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE: &GOTO -QUES : 
Stat 
fees) 30 FILES; 263 REC IN USE, 33 LEFT (of 296), 89% FULL (2 CYL) 
Mmeser * ftO2f002 
PLGe NOT FOUND 
ees. 13 LISTF * FTO2FO02 
men (00002) !!! 
Msc t * F*LOZFOOL 
PGeNAME FILETYPE MODE NO.REC. DATE 


SERG FTO2FOO1 Pl 53 9/16 
BERC ErUze Oo Ol Pi 5 9/16 
LNSERG Peo ZeOO Pi 3 9/16 
FILE PLOZEOO] Pil 3) 9/16 
ee ¢g f 


meeeos- NO RDR, NO PRT, NO PUN 
erase file ftO8f00l 

meeraine print file ftd02f001 
Betline print serg ftO2f00l 
offline print serc f£t02f001 
erase file frO2F001 

stat 

Sees, 26 FILES; 240 REC IN USE, 56 LEFT (of 296), 81% FULL (2 CYL) 
&goto -ques 

DO YOU WANT TO GO AGAIN? 

n 

BONTROL RETURNED TO CMS 

R; 


Figure 8. Sample user session with CMSWORK program 
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} 


called DATA FTO2F001, while the transformed series is 
written onto FILE FTO2F001. The parameters of the trans- 
form itself are written onto FILE FTO7FOOL. The TRANS 


program is self-contained, uSing no other programs. 


D. DIFF PROGRAM 

The DIFF program takes a given time series in FILE 
FTO2FO01, and performs non-seaSonal and/or seasonal differ- 
encing operations on the serieS in order to help achieve 
Stationarity for proper modeling. The user must indicate 
if seasonal differences are required; if so, the length of 
the season is input. The uSer is then required to enter 
the number of seasonal and nonseasonal differences. After 
Meecution of DIFF, the original data are in FILE FTO2FO001, 
and the differenced data are written onto FILE FTO3FO001. 
Geeuse the differenced data in any other program, CMSWORK 
program must be entered and the files altered to save the 
Original data (if desired) and to transfer the differenced 
data to FILE FTO2F001. The CMS commands to perform this 


memet ot files would look like: 


ALTER FILE FTO2FOOL Pl SAVEDATA FTO2FOOLl1 Pl 
R; 

ALTER FILE FTO3F00l1 Pl FILE FTO2FOOI1 Pl 

R; 


The DIFF program uses the IMSL subroutine called FTRDIF. 
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E. PLOT PROGRAM . 
The PLOT program plots any given time series which 
resides in FILE FTO2F001. Other than the time series, 
which must be on the user's disc or entered offline, the 
program requires only that an identification title for the 
plot itself be entered during execution. The plot is auto- 
matically printed offline. The PLOT program uses a modified 
PLOTP subroutine called PLOT8 from the IBM Scientific 


Subroutine Package Library (SSPLIB). 


F. AUTO PROGRAM 

The AUTO program takes a given time series in FILE 
FTO2F001 and calculates summary statistics of utility in 
Box-Jenkins modeling. The statistics include the sample 
mean, Variance, autocorrelations and partial autocorrela- 
tions for lags one though 40. Part of those statistics 
are printed out -at the user's terminal, while plots of the 
autocorrelations and partial autocorrelations versus lag 
are printed offline. These plots will also be written onto 
the user's disc in FILE FTO8F001, and may be plotted at 
the terminal in CMS outside the Editor or in the CMSWORK 
program; recall that thisis a time-consuming process. The 
AUTO program provides the user with essential information 
about stationarity, seasonality and model identification. 
With the Box-Jenkins procedures, the second moments (auto- 
correlations and partial autocorrelations) are valuable 


tools in tentative model identification. Additionally, AUTO 
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is useful in the diagnostic checkout phase of model 
building, when examining the model residuals for structure. 


The AUTO program uses the IMSL subroutine FTAUTO. 


G. ESTIMATE PROGRAM 

After the user has tentatively identified a model for 
a series through analysis of the series plot, autocorrela- 
tions and partial autocorrelations, ESTIMATE should be 
executed to calculate maximum likelihood estimates of the 
model parameters. ESTIMATE 1s written to calculate parameters 
for non-seasonal models only; for seasonal models, the 
programs YESTSEAS and WMARQRDT should be used for parameter 
estimation. The ESTIMATE program requires that the series 
to be modeled (already transformed and differenced to 
achieve stationarity) reside in FILE FTO2F001; additionally, 
the program will instruct the uSer to input the number of 
MA and AR parameters in the model, as well as the number of 
differences taken. The general model for which ESTIMATE 


calculates parameters is: 


o> (B) Ww, = 86 + 6(B) a, 


where we = voz! and Zy is the transformed value of Zh 
ESTIMATE calculates estimates of the AR parameters for the 
undifferenced form of the model, the estimated MA parameters 
and MA constant term, the residual variance, autocorrelations 


and partial autocorrelations of the residuals, plots of 
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these autocorrelations and partial autocorrelations, and 

a Chi-square goodness-of-fit statistic for the model. The 
original series to be modeled remains in FILE FTO2F001 

and the residuals of the model are stored in FILE FTO2F001. 
ESTIMATE uses the IMSL subroutine FTMAXL (a nonlinear 


gradient search algorithm) to estimate the model parameters. 


H. YESTSEAS PROGRAM 

The YESTSEAS program is used to calculate initial esti- 
mates of autoregressive and moving average parameters for 
a generalized seasonal Box-Jenkins model. These initial 
estimates are used as inputs to the program WMARQRDT, which 
then refines them to determine the parameter estimates that 
Minimize the sum of the squared residuals. The general 
model assumed in YESTSEAS is 


d 


o(B) 0(B°) 7°97 = 9(B)Q(B°)a 


De 
Se t 

The YESTSEAS program requires that the seasonal time series 
for which initial parameter estimates are required is located 
in FILE FTO2F001; it should be transformed, if desired, but 
generally not differenced, since YESTSEAS performs the 
required (seasonal and/or non-seasonal) differencing. 
Terminal inputs required by the program include the number 

of sdasonal and non-seasonal differences to be taken (if 


data is already differenced, enter zero for these values), 


the length of the seasonal period and the number of each 
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type of parameter to be estimated (that is, the number of 
each of the seasonal and non-seasonal AR and MA parameters). 
The output of the model, initial estimates for the requested 
parameters, is printed out at the terminal. YESTSEAS uses 
the IMSL subroutines FTRDIF, FTMAXL, FTAUTO and FTARPS, 


as well as a modified version of FTMAPS. 


I. WMARORDT PROGRAM 

The WMARQRDT program uses the initial parameter esti- 
mates calculated by YESTSEAS or other means as starting points 
for calculation of parameters of a generalized seasonal 
Box-Jenkins model which minimize the sum of the squared 
residuals. The general seasonal model assumed by WMARORDT 


can be written: 


$(B)O(B°)V°Voz, = @ + 6(B)O(B®)a,. 
Details of the Marquardt<-type non-linear least squares 
algorithm employed in this program are provided in Chapter 
III of this report. WMARQRDT requires that the series to be 
modeled reside in FILE FTO2F001; it should be transformed, 
if desired, but generally not yet differenced, since the 
program performs differencing. Inputs at the terminal 
include number of seasonal and non-seasonal differences, 
length of the seasonal period, number and type of parameters 
desired, and initial estimates of those parameters. The 


program then calculates the least-squares parameter estimates, 
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the standard error of these estimates, the MA constant term, 
the residual variance, the autocorrelations and partial 
autocorrelations of the residuals, the final model sum 

of squared residuals and a Chi-square goodness-of-fit sta- 
tistic (the portmanteau test). Additionally, it plots 
offline the autocorrelations and partial autocorrelations of 
the residuals. WMARQRDT uses the IMSL subroutines FTRDIF, 
FTAUTO, LINVIF and MDCDFI, the SSPLIB routine PLOTS8, and 

the Time Series Editor resident subroutines PARSH, MAROQORT, 
SUMSQ, SWAPB and FORMB. The calculations are performed in 


double-precision arithmetic. 


J. XSUMSQ PROGRAM 

The XSUMSQ program accpets an already transformed and 
differenced series in FILE FTO2F001 and general seasonal 
and non-seaonal Box-Jenkins parameter values for a general 
seasonal model from the terminal. The program then calcu- 
lates the residual sum of squares for these parameters, 
giving the user a feel for the "goodness-of-fit" the 
parameter estimates provide in modeling the given time 
series. XSUMSQ uses the Time Series Editor resident sub- 
program XSUMSQ. All output is printed at the terminal. 


The calculations are performed in double-precision arithmetic. 


K. FORECAST PROGRAM 
The FORECAST program uses the estimated Box-Jenkins 


(Seasonal or non-seasonal) model to compute forecasts of 





the transformed and differenced time series. FORECAST re- 
quires that the time series to be forecasted reside in FILE 
FTO2FOO1l in transformed form. During execution the input 
required at the terminal includes the forecast origin, the 
numbers and estimated values of the AR and MA parameters 

of the model, the overall MA constant, the plot origin 
index, the maximum forecast lead time, the order of differ- 
encing in the model, and a level of mig nmrcance for the 
forecast confidence limits. The program needs the transfor- 
mation parameters from FILE FTO7FO01, created when using 
program TRANS, and performs the inverse transform to return 
the data to its original form. FORECAST also undiffer- 
ences differenced data. The program output includes the 
forecasts up to the given maximum lead time, the deviations 
from each forecast for the confidence limits, and plots of 
forecasts and confidence limits plotted offline. The 
FORECAST program uses IMSL subroutine FTCAST and the special IBM 
pOPLIB Subroutine UTPLT8. The plot itself and output data 
are written onto FILE FTO8F0O01 on the user's disc; since 
this plot file takes considerable space, it is usually 

best to erase it after execution of the program is completed 


and the plot itself has been printed offline. 


L. ROOTS PROGRAM 
The ROOTS program is used to determine the roots of 
the characteristic equation for a general ARIMA model. The 


Beers are useful for testing for stationarity and for 
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determining the gorau of the forecast function. ROOTS uses 
the IMSL subroutine ZPOLR to calculate the roots; it requires 
as input the number of AR parameters in undifferenced form,~ 
and the actual values of these parameters, as determined 

from the ESTIMATE or WMARQRDT programs. The roots are 
printed out at the terminal. If desired, roots can also 

be obtained for the MA polynomial, simply by substituting 


the appropriate MA values for the AR values as input. 


M. GENERATE AND SIMULATE PROGRAMS 

Because of their similarities, the GENERATE and SIMULATE 
programs are described together. The GENERATE program per- 
mits the user to generate a time series from any non-seasonal 
ARIMA model he specifies. The user must identify the model 
and give values for the model parameters and starting con- 
ditions; a random number seed must also be input. The 
program takes the specified model, generates random noise 
terms, and calculates as many values of the time series as 
desired. The GENERATE program can be useful for purposes 
ofd classroom instruction in the generation of a wide variety 


of time series examples for model identification exercises. 


tcuppose the model was ARIMA(1,1,0). The differenced form 


would be Bon) Vz) = ae -- a,: Then, the undifferenced form 
2 


would be written (1-(1+9,)B + >,8 )2, = oe + ay , Where the 
ul 


left-hand side was found by multiplying (1-$,8) By eel B) aS, 
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It could also be useful in the diagnostic phase of model 
checkout; a time series could be generated from the esti- 
mated model, and itS properties compared with those of 

the original series. If large discrepancies occur in this 
comparison, this may be evidence that the model is inade- 
quate. The model output, the generated time series, is 
written in FILE FTO02F001 as well as being printed offline. 
The GENERATE program uses IMSL subroutine FTGENI1. 

The SIMULATE program provides the capability of gen- 
erating any number of simulated time series. This program 
is useful for prediction of what might happen in the future; 
it also demonstrates that within a given model, the time 
series actually observed can vary considerably. This program 
uses GENERATE, but also requires as input the actual series 


1 and the number of simulated series the 


in FILE FTO2F001 
user wishes to generate. The series should be already 
transformed and differenced. Additionally, SIMULATE allows 
the user to select values of the original time series as 

| Starting values for the simulated series. The output con- 


Sists of the simulated series printed at the terminal. The 


SIMULATE program uses IMSL subroutine FTGENI1. 


me «=ORELP OPTION 
The HELP option is a modification to the Time Series 


Editor CP/CMS program itself that allows the user to obtain 


lie desired, the user may choose to input his own starting 
values, and not use the actual time series values. 
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information about Editor programs after the initial 


introductory phase has been completed. 


O. SUMMARY 

This chapter has provided brief descriptions of the 
options available for analysis of time series data in the 
Time Series Editor. The following two chapters contain 
detailed examples of time series analysis for both a non- 
seasonal and a seasonal series; the accompanying terminal 
sessions and computer output for analysis of these series 
are contained in full in Appendices B and C. Appendix A 
to this report contains a User's Guide for the Tee ce cies 
Editor, while Appendix D contains a complete listing of 
the Time Series Editor and resident programs. Appendix E 
contains listings of the data sets used in the ara ecs of 


the next two chapters. 
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V. EXAMPLE NON-SEASONAL TIME SERIES ANALYSIS 


This chapter provides a description of a complete analy- 
sis of a non-Seasonal time series uSing the Time Series 
Editor. The time series analyzed is Series C (Chemical 
Process Temperature Readings), taken from Box and Jenkins 
[Ref. 4, p. 528]. This time series was selected since it 
is fully analyzed and discussed in -Reference 4, so that the 
interested user can compare for himself the results pre- 
sented there with the results given by the Time Series 
Editor. After reading this overview of the analysis, the 
reader can follow the actual user session and associated 
user session Editor program output contained in Appendix B 
Bomtnis report. 

The first step in the analysis of series C 1s to plot 
the series uSing program PLOT. The plot reveals rather wide 
fluctuations in the series, but not the sort of explosive 
nonstationary behavior that would render an attempt to model 
Beuecress. The plot also reveals that the time series has 
a large amount of momentum (movements of the series tend to 
resist changes of direction over time). This is character- 
istic of ARIMA models with two or more differences. 

As a second step, the autocorrelations (autos) and 
partial autocorrelations (pautos), mean, and variance of 
the series were estimated uSing AUTO. The plot shows that 


the autos dampen out slowly in a nearly linear fashion. 
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This is an indication that the process is nonstationary 
and that one or more differences are required to make it 
stationary. The pauto plot is not informative when the 
autos fail to dampen out rapidly. 

As suggested by the plots of the original series and 
its auto and pautos, the program DIFF was used to examine 
the series with both one and two non-seasonal differences; 
i.e., the series zy = Viz, and Zy = Vez, were derived. 

The plots suggest that two differences might be required to 
achieve stationarity in the series. However, one might be 

able to get by with a single difference and an AR parameter 
near unity, since the autos of z!' dampen out slowly. Thus, 


c 


two possible model candidates are suggested: 


MeArrMA(1.1,0): (1 - > B)V 2, = 9 +a, , 


2 = 2 
free ARIMA(0,2,2): V Zz, = Oe (1 8,8 9.8 a, 
For purposes of estimation, the second model was extended 
to two moving average parameters. Such overfitting is often 
done to see if the estimated. moving average parameters turn 
Cut to be near zero, thus confirming the tentative identification. 


The next step in the analysis process is to calculate 


the maximum likelihood estimates of the model parameters. 
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The program ESTIMATE was used to do this. The estimated 


parameters for the ARIMA(1,1,0) model were: 


6, = 0.8073 
8 9 = -~0,006789 
a* = 0.0177633 


The autos and pautos of the residuals were calculated 
to test the model. The correlations should appear to be 
estimates of a pure white noise process if the model is 
adequate; the (1,1,0) model seemed to pass that test. The 
chi-square goodness-of-fit test results ina ve ee or Go cmawele m) 
24 df and a signficance level of 0.2247. Thus, there is 
no strong evidence tec suggest that the (1,1,0) model is 
inadequate. The model can be written, in undifferenced 
orm : 


(1 - 1.8073B + 0.8073B°)z, = a, . 


The eS 1s omitted here since it is nearly zero in value. 


The estimated parameters for the (0,2,2) model were 


aT Seo s373 
05 = 0.1296 
- = -0.0026787 
of = 0.0189559 


As before, the correlation plots of the residuals fail to 
Suggest any inadequacy of the model. However, the chi-square 
goodness-of-fit test for this model yielded a ve = 36:75 


with 23 df and a significance level of 0.0346. Therefore, 
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if the ARIMA(0,2,2) model were a "correct" one, a chi-square 
value as large as 36.75 would occur by chance with a proba- 
ity of only 0.0346. Due to its simplicity (parsimony 

is a desirable feature in time series models) and better 
me, the ARIMA (1,1,0) model was chosen as the "“best" of 

the two alternatives originally selected. 

The ARIMA(1,1,0) model calculated was then used to make 
forecasts and confidence limits for these forecasts, using 
program FORECAST. The forecasts were made for 25 periods 
into the future, with plot origin 200 and forecast origin 
220. The length of the series C itself is 226; with one 
difference, it becomes 225. The plot in Appendix B shows 
the forecasted values, along with the confidence limits at 
level of significance 0.10; the plot shows that the width 
of the confidence limits increases substantially as the lag 
gets large. 

Finally the program ROOTS was used to calculate the 
roots of the characteristic equation for the ARIMA(1,1,0) 
model, which can be written: 


ie 073° = 0.8073B- = 0 


The roots calculated by ROOTS were 1.239 and 1.0. Also, 

the program SIMULATE was used to generate a simulated series 
from the ARIMA(1,1,0) model, using the last two values of 
the actual series as starting values. The entire session 
lasted about two hours, including checkout of plots and 


consumed less than a minute of CPU time. 
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VI. EXAMPLE SEASONAL TIME SERIES ANALYSIS 


This chapter provides a description of a complete analy- 
Sis of a seasonal time series using the Box-Jenkins tech- 
nique through the Time Series Editor. The time series 
analyzed is Series G (International Airline Passengers: 
Monthly Totals (Thousands of Passengers) January 1949 - 
December 1960), taken from Box and Jenkins [Ref. 4, p. 531]. 
This time series was selected since it was fully analyzed 
and discussed in Chapter 9 of Reference 4; again, the 
interested reader can compare for himself the results pre- 
sented by Box and Jenkins with those given by the Time 
Series Editor programs. After reading this overview of the 
analysis, the reader can follow the actual user session and 
associated Editor program output contained in Appendix C 
Boechnils report. 

The plot shows both annual seasonality and an increasing 
trend. This suggests that a seasonal model with differencing 
May be needed. The autocorrelations (autos) and partial 
autocorrelations (pautos) were then calculated and plotted 
using AUTO. The plot of the autos shows very slow damping 
Out with peaks at intervals of 12, indicating strong non- 
stationarity and a seasonal period of 12 months. Again, 
here the pauto plot is not informative, due to the slow 
damping out of the autos. As it is sometimes useful in 


dealing with seasonal models, a natural log transform was 
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made on the series G data uSing program TRANS. The plot 
of the logged data shows little change in non-stationarity, 
and the autos and pautos of the logged data still indicate 
seasonality and non-stationarity. 

Next, differencing to achieve stationarity was tried, 
using program DIFF; both one seasonal and one non-seasonal 
difference of the logged data was taken. The autos and 
pautos of the data WV 5 


rapidly damping out. They also exhibit sharp peaks at the 


(ln Z,) appear much improved, 


periods of seasonality, indicating the need for seasonal 
parameterization in the model. Since positive identification 
of the model type from the auto and pauto plots was not 


possible, two were postulated as close candidates: 


& 


Ih 


meee ARLMA (0,1,1) x (0,1,1) a weleten: 

Pes = _ = 12 

V Vy (in Z,) = 8G + (l 9,8) (1 ©, 128 a, P 

and 
12 ' 
Yee ARIMA ({1,1,1) x (0,1,1) , written: 
= 151 x a 7 1,2 
fl: 8,B)V Vio (in Z,) = - + (1 9,8) (1 oF ie )a,- 


The next step in the analysis was to calculate the 


parameter estimates for both models. YESTSEAS was used to 
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obtain initial parameter estimates. For the model 


ARIMA Cee (onl 1) --, YESTSEAS yielded starting values: 


Nn 


ue) 0.390425 


Te ais) = 0.534397 . 
For this model, calculation of final parameter estimates 
in WMARQRDT using the starting values above yielded the 


following in seven iterations: 


ores One 7 7152 
Sag = OoSvaEeu 
6 = 0.000291 ; 
© 


the chi-square statistic for residual lack of fit of the 
model was a = 29.713571, with df = 38 and a probability 
of exceeding the ie value of .829. 

For the ARIMA (1,1,1) x (0,1,1)+* model YESTSEAS calcu- 


lated starting values: 


PL hss Helizeg9 


AN 


0 0.490883 


1 8) 
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91 ,12(s) Sig Bie) 219 7/5) 


for this model, WMARQRDT yielded the following parameter 


Palues in ten iterations: 


“w 


o, = 0.146667 
84 = 0.509816 
2 ag = WeBussels 


The chi-square statistic for residual lack of fit of the 
model was ve = 28.302543, with df = 37 and a probability 
of exceeding the ve value of .847. 

Next, the plots of autos and pautos of the residual for 
both models were examined; both models have autos and pautos 
that dampen out rapidly, and essentially appear to be 
"white noise", with the exception of a small peak at auto 
and pauto of lag 12, corresponding to the seasonal period. 
The chi-square statistics indicate that both models have 
"about" the same goodness of fit to the data. Consequently, 
both models were used to make forecasts of the series. 

Program FORECAST was used to forecast series G into the 
future for 30 time periods, with a forecast origin of 131; 
the plot origin selected was index 100. As expected, both 


models do about the same job of forecasting the series. 


ou 





B@eesignifticance level for the confidence limits for both 
plots was 90%. 

As a final check of both models, the program XSUMSQ 
was used to calculate the sum of squared residuals for 


models. The results were as follows: 


ae For ARIMA (On 1,.1) x Gene) r 


Nn a2? r 
eZ 
meer or ARIMA (1,1,1) x (0,1,1) ; 


nan wN wN 


“A 2 - 
S (dy ee 1! = 0.181643 


It can be seen here that both models are nearly equally 
"good", and either could be used with a reasonable amount 


of confidence for forecasting series G. 
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VII. SUMMARY 


There has been an increasing need in military applica-. 
tions, aS well as in the general business area, for a 
computer-aided time series analysis and forecasting capability 
that is both powerful and easy to use. Current military- 
associated functions in which time series techniques are 
valuable include material inventory Management, recruiting, 
personnel management, budget analysis, force level projections, 
and perhaps even short-term weather prediction. 

The Time Series Editor that has been described in this 
report provides a unified collection of programs that greatly 
facilitate the conduct of a complete time series analysis 
and forecasting evolution using Box-Jenkins methodology. 

The routines guide the user from the data input stage, 
through the data analysis and parameter estimation stages, 
and finally through the model building, diagnostic checking 
and forecasting stages. With its simple input requirements, 
the Editor can be quickly mastered by even a beginning 
computer user, who has a basic understanding of Box-Jenkins 
methodology. 

The Box-Jenkins methodology has been described broadly 
in Chapter II and descriptions of the programs resident 
in the Time Series Editor have been given in Chapter IV. 
Example analyses of both non-seasonal and seasonal time 


series are provided in Chapters V and VI, with complete 
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user sessions and output packages for these analyses 
included as Appendices B and C. The Appendices also con- 
tain a Time Series Editor User's Guide, as well as complete 
program listings. 

Although the Time Series Editor in its present form 
provides an excellent capability for interactive time 
series analysis uSing the Box-Jenkins methodology, there 
are additions to the program set that might be made in the 
future which would improve the Editor's capability and 
utility. It is first recommended that the model diagnostic 
testing capability be extended to include a periodogram 
analysis, and/or other tests related to the spectral analy- 
Sis of time series {see Ref. 4, p. 294]. Another useful 
addition might be to include in the Time Series Editor a 
package capable of improving the ease of entering data 


and simplifying the input error correction problem. 
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APPENDIX A 


USER'S GUIDE TO THE TIME SERIES EDITOR 


The Time Series Editor is a collection of FORTRAN 
programs driven by a control program called TIMESER EXEC 
in the CP/CMS executive language that has been designed 
specifically for the analysis and forecasting of time 
series data using the general Box-Jenkins methodology. 
This Guide contains information essential for the user to 
access the Time Series Editor, enter time series data, 


build a model and obtain the desired output. 


I. DATA INPUT 

The Time Series Editor requires that time series data 
be entered into the sequential FORTRAN input/output file 
named FILE FTO2FOO1l. This can be done either online or 
via cards read offline. The following diagram illustrates 
the proper card deck arrangement. When your deck is ready 
for input, give it to a computer center system operator 
memdirect it to your disc space. 

If the user has a data deck already punched up ina 
format other than 5F15.6, he may enter the series as above 
into FILE FTO3F001 (without the length of series card) and 
use the Editor program ZFORMAT to transform it into the 
FILE FTO2FOO01 in the proper format, without destroying his 


Original file. This program is described in the next section. 
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i:. TABLE OF OPTIONS 
This Table provides the user with the basic information 
necessary to understand the data requirements, functions 


and output options for each program in the Editor. 
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=e. THE BASIC USER SESSION 

To use the Time Series Editor, the user must log into 
CMS, get into CP, link to the disc storage area where the 
Time Series Editor resides, reimplement CMS, log into the 
general user and Time Series Editor disc areas, and enter 
the TIMESER routine. This section will provide explicit 
guidelines to enable the user to perform the above steps 
on the NPS CP/CMS system. Commands marked with an asterisk 
(*) are those actually entered on the terminal by the user 
(the asterisk itself is omitted). Those without an aster- 
isk and those written in all capital letters are system 
responses at the terminal. Numbered sentences are comments, 
which will not appear during an actual user session. The 
instructions and system responses assure the user is on an 
IBM 2741 Input/Output Terminal. Some minor modifications 


may be necessary if other terminals are used. 


1. Turn the terminal on, depress the BREAK key, and wait 
for the system to respond: 


CP-67 online xd.65 gsyosu 


2. Depress the ATTN key. The roll bar will advance and 
the keyboard will unlock. Then enter: 


*login aaaapbb 450k 


3. aaaa is the user's identification number, and nn is the 
terminal number (usually written on the terminal). 
For example, if the user's ID number is 1621 and the 
terminal number is 44, the input would be: 
Meqgin 1621p44 450k. The addition of 450 k to the normal 
login command is necessary to execute the program 
WMARORDT in the Editor; for users not planning to 
execute this program during a session, this addition 
is not necessary. 


LO5 





7G . 


#1. 


#2. 


i. 


14. 


The system will respond with the statement: 
ENTER PASSWORD: 


The user then enters his password, or the general users 
password npg; 


*password 
The system will then respond: 


ENTER 4-DIGIT PROJECT NUMBER FOLLOWED BY 4-CHARACTER 
Seal CENTER CODE: 


The user then enters: 
eaggghhhh 


gggg is the asSigned project number, and hhhh is the 
user's Section deSignator or the faculty code. 


The system will respond with the message of the day, 
Such as: 


Mees HOURS ... 0930=2200 (MON-THURS) ... 0930-1800 (FRI) 


‘OUTPUT RETAINED 5 DAYS 


Cms Version 3.25 

t this point, the user 1s in CMS. He must then get into 
CP; this can be done by hitting the ATTN key. The system 
will then respond: 

CP 


The user must then link to the TIME SERIES EDITOR; 
this 1S accomplished by entering: 


wank 2069p 191 192 

The system will respond with: 

ENTER PASSWORD: 

The password (read only) to enter the Editor is: 
*timser 

The system then responds: 


SET TO READ ONLY 
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15. 


16. 


m/. 


18. 


19. 


20. 


zL. 


ae 


Zz . 


The user now implements CMS by: 
*ipl cms 

The system will respond: 

CMS Version 3.25 


Now the user must log into both the general user and 
the Time Series Editor area by entering: 


Blogin 191 

The system vpialal respond with a message such as: 
R; 

The user then enters the command: 

miogin 192 t,p 


The system will respond: 


me (192) R/O 


R; 


The user can then enter the Time Series Editor 
(guided version) by entering the command: 


*timeser 

The system will respond: 

EACH 2 SECONDS EXECUTION TIME IS INDICATED BY * 
OU HAVE ENTERED THE TIME SERIES EDITOR 

PLEASE RESPOND TO EACH QUERY WITH AN INPUT AT THE TERMINAL. 
ENTER ONLY THE FIRST LETTER FOR A WORD RESPONSE. 

ENTER NUMERICAL VALUES VIA FORTRAN FORMAT. 

TYPE INTEGER VALUES (RIGHT JUSTIFIED) FOR NAMES STARTING 
WITH I THROUGH N. TYPE FLOATING VALUES WITH DECIMAL FOR 
ALL OTHERS. 

DO YOU WANT A LIST OF THE OPTIONS? 

The user is then on his own, guided by the Exec routine. 
See the notes that appear at the end of this guide for 
additional information. Eventually the user will be 


asked: 


B@,yYOU WANT TO TRY AGAIN? 
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24. If a yes response is given, another sequence will begin; 
if the response is no, the user will be taken out of 
the Time Series Editor environment and returned to CMS. 
The system response will be: 


CONTROL RETURNED TO CMS 
R; 


mee) 6The: «user can then log out of CMS by typing: 
eep logout 
26. The system will respond with: 


CONNECT= 00:08:02 VIRTCPU= 000:07.98 TOTCPU= 000.10.94 
LOGOUT AT 14.22.04 on 10/16/78 


27. The user should then turn off his terminal and tear off 
the output from his session. 


The more experienced user can dispense with the 
“welcome aboard" section of the Time Series Editor and get 
right down to business by using the shortened version of the 
Editor. This shortened version may be entered by linking in 
the normal way, and then entering the Editor by typing the 
COMMAND 
*timeser s (asterisk omitted). 
The system will immediately respond: 


ENTER LETTER FOR OPTION YOU WANT. 


The session inside the Editor then begins. 
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IV. BRIEF SAMPLE USER SESSION 


A brief sample user session is given below; it includes 


copies of the offline output generated during the session. 


G2 


space 





rec, a beef ile = 


ak a 





repeat login nur@pné@ 


login 1621p44 450k 
ENTER PASSWORD: 


ENTER 4-DIGIT, PROJECT NUMBER FOLLOWED BY 4-CHARACTER COST CENTER CODE: 
0444r172 | s 

READY AT 17.24.38 ON 09/16/78 

CMS Version 3.25 


Stat 

P (191): 29 FILES; 241 REC IN USE, 55 LEFT (of 296), 81% FULL (2 CYL) 
R} 

eo gf 

FILES:~ NO RDR, NO PRT, NO PUN 

R; 

CP 

link 2069p 191 192 

ENTER PASSWORD: 


SET TO READ ONLY 


ipl cms 
CMS Version 3.25 


login 191 


R; 

login 192 t,p 
T (192) R/O 
R; 


timeser 
BACH 2 SECONDS EXECUTION TIME IS INDICATED BY *. 


YOU HAVE ENTERED THE TIME SERIES EDITOR. 

PLEASE RESPOND TO EACH QUERY WITH AN INPUT AT THE TERMINAL. 
ENTeR ONLY THE FIRST LeTTER FOR A WORD RESPONSE. 

ENTER NUMERICAL VALUES VIA FORTRAN FORMAT. 


TYPE INTEGER VALUES (RIGHT JUSTIFIED) FOR NAMES STARTING 
WITH I THRU N. TYPE FLOATING VALUES WITH DECIMAL FOR ALL OTHERS. 


DO YOU WANT A LIST OF THE OPTIONS? 


Y 

OPTION DESCRIPTION 

GENERATE ----- GENERATE ANY ARIMA TIME SERIES 

AUTO --------- CALCULATE AUTOCORRELATIONS, PAUTOS, MEAN AND VARIANCE 
PLOT ---------PLOT A TIME SERIES 

ESTIMATE ----- CALCULATE MAX LIKELIHOOD ESTIMATES OF ARMA PARAMETERS 
DIFF --------- DIFFERENCE A TIME SERIES 

FORECAST ----- FORECAST FUTURE VALUES, CONSTRUCT CONFIDENCE INTERVALS 
TRANS -~------ “TRANSFORMS VALUES OF A TIME SERIES 

ROOTS -------- DETERMINES ROOTS OF ARIMA CHARACTERISTIC EQUATION 
ZFORMAT ------ ALTER DATA FILE TO FORMAT 5F15.6 

CMSWORK ------ PERFORM CP/CMS COMMANDS IN TIMESER EXEC 

SIMULATE -----SIMULATE NONSEASONAL TIME SERIES 

YESTSEAS -----CALCULATE INITIAL SEASONAL PARAMETERS 

@M4ARQRDT ----- MARQUARDT SOLUTION FOR PARAMETER ESTIMATES 

XSUMSQ ------- CALCULATE SUM OF SQUARES FOR ARBITRARY PARAMETERS 


WOULD YOU LIKE MORE INFO? 
“4 


1A 





ENTER OPTION YOU WANT INFO ABOUT. 
a 


oOo THIS PROGRAM CALCULATES AUTOCORRELATIONS, PARTIAL 
AUTOCORRELATIONS, THE MEAN AND THE VARIANCE FOR A GIVEN TIME SERIES 
wHICH MUST RESIDE IN FILE FTO2F001. THE PROGRAM USES 

FTAUTO IN THE IMSL LIBRARY. THE AUTOCORRELATIONS AND PAUTOS CAN BE 
PLOTTED OFFLINE. 


DO YOU WANT INFO ABOUT ANOTHER OPTION? 


n 

DO YOU WANT TO TRY A SESSION? 

Y 

ENTER LETTER FOR OPTION YOU WANT. 
Cc 


ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE: &GOTO ~-QUES 

Mmiste * £t02f001 

FILENAME FILETYPE MODE NO.REC. DATE 


SERG FTO2FOOl Pl 3 9/16 
SERC FTO2FOOl Pl 5 9/16 
LNSERG FTO2FOO1l Pl 3 9/16 
FILE BEOZFOO! Pl 1 9/16 


erase file f£tO2f£001 

alter sere f£t02f£001 pl file f£t02£001 pl 

stat 

P (191): 28 FILES; 240 REC IN USE, 56 LEFT (of 296), 81% FULL (2 -CYL) 
&goto ~ques 

DO YOU WANT TO GO AGAIN? 

y 


ENTER LETTER FOR OPTION YOU WANT. 
a 
IS YOUR DATA IN FILE FTO2F001? 


Y 

EXECUTION BEGINS... 
AUTOCORRELATIONS 
e 


0.978 0.944 0.902 0.854 0.802 0.748 0.692 0.635 0.579 0.323 
0.468 0.413 0.359 0.305 0.253 0.201 0.150 0.098 0.047 -0.003 
205052 -0.101 -0.15] -0.200 -0.248 


PARTIAL AUTOCORRELATIONS 


0.978 -0.260 -0.157 -0.093 -0.058 -0.045 -0.012 -0.038 -0.022 ~-0.010 
-0.036 -0.041 -0.038 -0.024 -0.037 -0.027 -0.032 -0.070 -0.048 -0.024 
=0.034 -0.061 -0.079 -0.048 -0.037 


MEAN= 22.9739 VARIANCE = 4.22273 


ENTER TITLE FOR PLOTS. 

autos and pautos for series c data 

YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 

DO YOU WANT TO GO AGAIN? 


Y 
ENTER LETTER FOR OPTION YOU WANT. 


o4 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE: SGOTO -QUES 

offline print file ft0O2£001 

alter file £t0O2f001 pl sere fFtO2f001 pl 
&GOtoO ~ques 

DO YOU WANT TO GO AGAIN? 


n 
CONTROL RETURNED TO CMS 
R; 
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V. PROBLEM CONTROL NOTES 

This section will cover corrective measures that can be 
taken when things fail to go as expected while in the 
TIMESER environment. 

a. A typing error in the CMS environment can be 
corrected by typing the @ character as many times as is 
required to back up and then type the correct values. For 
example, if the user typed timesre, the user could correct 
the mistake by typing two @ signs, followed by the correct 
spelling er, as follows: timesre@@re. An entire line can 
be deleted by typing the ¢ character (or [ On some terminals). 

b. When working with TIMESER executive programs, the 
user should exercise care before hitting the return key (or 
control Ss on some terminals). If an input value is required 
and the return key is hit before the proper response is 
entered, the user will be likely to get thrown out of the 
editor and have to begin again. In most cases, errors can 
be corrected only before the return key is struck (in some 
programs you get a second chance for input). 

c. Particular care should be taken for integer value 
input, which must be right justified in the format field. 
The editor will advise the user in all cases where the 
integer format is other than Il. 

d. If for any reason the user finds himself in a debug 
Or error condition (caused by erroneous data, or a "blowup" 


in one of the non-linear optimization routines usually 


1S 


caused by very poor initial input values), the following 
procedure will get the user back into the normal TIMESER 
environment: 
(1) depress the ATTN key twice; this gets 
the user into CP; hit the ATTN key again, 
and then type kx to kill the execution. 
(2) re-ipl CMS, and login 191 and then login 
Oza, D - _— 
(3) then the user can type TIMESER or TIMESER S, 
and return to the TIMESER environment 
On the next page 1S a sample user session where an error 
cauSing a debug condition has occurred. 
Here the user executed a program that required data 
Pee LLE FTO2FO0O1, and the FILE did not exist. As the 
example shows, recovery iS quick. Simply hit the break 


econ tO get into CP, login 191, login 192 t,p, and then 


type timeser s to return immediately to the Editor environment. 
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APPENDIX B 
KEYBOARD SESSION AND OUTPUT FROM ANALYSIS OF 
TIME SERIES C (NON-SEASONAL) 

This appendix contains the complete terminal keyboard 
session and all offline output from the analysis of Box 
and Jenkins' [Ref. 4] time series C (Chemical process 
temperature readings) uSing ning, aelnte Series Eater. “This 
appendix is intended to supplement the discussion of the 
analysis of series C discussed in Chapter V of this report. 

The terminal listings are presented first, followed by 


the computer offline output. 


ie? 





repeat login nuré@pn@ 


login 1621p44 400k 
ENTER PASSWORD: 


ENTER 4-DIGIT PROJECT NUMBER FOLLOWED BY 4-CHARACTER COST CENTER CODE: 
0444r172 

READY AT 14.23.02 ON 09/16/78 

CMS Versior 3.25 


Piste * £t02f001 
FILENAME FILETYPE MODE NO.REC. DATE 


FILE FTO2ZFOO1L Pl 3 9/15 
SERG FTO2FOOl pl 3 9/16 
SERC FTO2FOO1l pl 5 9/16 

R; 

erase file ft02£f001 

R; 

alter serc f£t02f001 pl file £tO02£001 pl 
R; : 


CP 
link 2069p 191 192 
ENTER PASSWORD: 


SET TO READ ONLY © 


ipl cms 
CMS Version 3.25 


login 191 
R; 
login 192 t,p 
T (192) R/O 
R; 

wer 
~ACH 2 SECONDS EXECUTION TIME IS INDICATED BY *. 


YOU HAVE ENTERED THE TIME SERIES EDITOR. 

PLEASE RESPOND TO EACH QUERY WITH AN INPUT AT THE TERMINAL. 
ENTER ONLY THE FIRST LETTER FOR A WORD RESPONSE. 

ENTER NUMERICAL VALUES VIA FORTRAN FORMAT. 


TYPE INTEGER VALUES (RIGHT JUSTIFIED) FOR NAMES STARTING 
WITH I THRU N. TYPE FLOATING VALUES WITH DECIMAL FOR ALL OTHERS. 


DO YOU WANT A LIST OF THE OPTIONS? 


Y 

OPTION DESCRIPTION 

GENERATE ----- GENERATE ANY ARIMA TIME SERIES 

AUTO ---e<---- CALCULATE AUTOCORRELATIONS, PAUTOS, MEAN AND VARIANCE 
PLOT w-2------ PLOT A TIME SERIES 

Zc LMATE -<---- CALCULATE MAX LIKELIHOOD ESTIMATES OF ARMA PARAMETERS 
DIFF s--32e----- DIFFERENCE A TIME SERIES 

FORECAST <---- FORECAST FUTURE VALUES, CONSTRUCT CONFIDENCE INTERVALS 
TRANS wo3rrre--- TRANSFORMS VALUES OF A TIME SERIES 

es ———————— DETERMINES ROOTS OF ARIMA CHARACTERISTIC EQUATION 
ZFORMAT <------ ALTER DATA FILE TO FORMAT 5F15.6 

CMSWORK o----- PERFORM CP/CMS COMMANDS IN TIMESER EXEC 

SIMULATE -----SIMULATE NONSEASONAL TIME SERIES 

YESTSEAS ----- CALCULATE INITIAL SEASONAL PARAMETERS 

WMARQRDT ----- MARQUARDT SOLUTION FOR PARAMETER ESTIMATES 

XSUMSQ w------ CALCULATE SUM OF SQUARES FOR ARBITRARY PARAMETERS 


WOULD YOU LIKE MORE INFO? 
Y 


P20 





ENTER OPTION YOU WANT INFO ABOUT. 
a 


ao) -—-<----— THIS PROGRAM CALCULATES AUTOCORRELATIONS, PARTIAL 
AUTOCORRELATIONS, THE MEAN AND THE VARIANCE FOR A GIVEN TIME SERIES 
WHICH MUST RESIDE IN FILE FTO2F001. THE PROGRAM USES 

FTAUTO IN THE IMSL LIBRARY. THE AUTOCORRELATIONS AND PAUTOS CAN BE 
PLOTTED OFFLINE. 


DO YOU WANT INFO ABOUT ANOTHER OPTION? 


y 

ENTER OPTION YOU WANT INFO ABOUT. 

d 

a oo THIS PROGRAM TAKES A SERIES IN FILE FTO2F001, ASKS IF 


IT IS SEASONAL, AND THEN PERFORMS NONSEASONAL AND/OR SEASONAL 
DIFFERENCING, AS REQUESTED. THE ORIGINAL DATA REMAINS IN FILE 
FTO2FO001, WHILE THE DIFFERENCED SERIES IS PUT IN PILE FTO3FOO1. 


é 


DO YOU WANT INFO ABOUT ANOTHER OPTION? 


n 
DO YOU WANT TO TRY A SESSION? 
Y 


ENTER LETTER FOR OPTION YOU WANT. 
a 
IS YOUR DATA IN FILE FTO2F001? 


Y 
EXECUTION BEGINS... 
AUTOCORRELATIONS 


0.978 0.944 0.902 0.854 0.802 0.748 0.692 0.635 0.579 0.523 
Meroe, 0.413 0.359 0.305 0.253 0.201 90.150 0.098 0.047 -0.003 


=0.052 -0.101 -0.151 -0.200 -0.248 


PARTIAL AUTOCORRELATIONS 


0.978 -0.260 Boe si -0.093 -0.058 -0.045 -0.012 -0.038 -0.022 -0.0)30 


-0.036 -0.041 -0.038 -0.024 -0.037 -0.027 -0.032 -0.070 -0.048 -0.024 


-0.034 -0.061 ~0.079 -0.048 -0.037 


MEAN= 22.9739 VARIANCE = 4.22273 


ENTER TITLE FOR PLOTS. 

autos and pautos of series c data / undifferenced 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM I140 UNDER YOUR USER ID NUMBER. 

*DO YOU WANT TO GO AGAIN? 


y 


Idk 





ENTER LETTER FOR OPTION YOU WANT. 
a : : 
IS YOUR DATA IN FILE FTO2F001? 


Y 
EXECUTION BEGINS... 
IS YOUR TIME SERIES SEASONAL? 


n 

ENTER NUMBER OF NONSEASONAL DIFFERENCES. 

1 

*DO YOU WANT TO PLOT AUTO AND PAUTO OF TRANSFORMED DATA? 


y 
EXECUTION BEGINS... 


AUTOCORRELATIONS 

0.805 0.653 0.526 0.442 0.380 0.318 0.262 0.186 0.139 
0.097 0.094 0.074 0.073 0.070 0.072 0.089 0.048 0.041 
6.044 0.048 0.000 -0.052 -0.097 

PARTIAL AUTOCORRELATIONS 
0.805 0.010 -0.007 0.051 0.027 -0.019 -0.013 -0.080 0.020 
~-0.137 @.094 -0.027 0.034 0.006 0.012 0.043 -0.121 0.060 
0.038 


-0.022 -0.116 -0.076 -0.008 


MEAN=-.346667E-01 VARIANCE = 0.531982E-01 


ENTER TITLE FOR PLOTS. 

autos and pautos of series c data / one nonseasonal difference 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 

PICK UP IN ROOM I140 UNDER YOUR USER ID NUMBER. 


DO YOU WANT TO GO AGAIN? 
y 


ENTER LETTER FOR OPTION YOU WANT. 
d 
IS YOUR DATA IN FILE FTrO2F001? > 


y 
EXECUTION BEGINS... 
IS YOUR TIME SERIES SEASONAL? 


n 

ENTER NUMBER OF NONSEASONAL DIFFERENCES. 

Z 

DO YOU WANT TO PLOT AUTO AND PAUTO OF TRANSFORMED DATA? 


Y 
EXECUTION BEGINS... 


AUTOCORRELATIONS 
-0.079 -0.065 -0.122 -0.063 0.013 -0.018 0.049 =-0.052 =-0.124 
“0.122 0.072 -0.077 0.029 -0.011 -0.058 0.171 -0.101 -0.013 
-0.007 0.135 0.014 -0.013 -0.135 
PARTIAL AUTOCORRELATIONS 

-0.079 -0.072 -0.135 -0.094 -0.022 -0.051 0.022 -0.060 -0.145 
-0.143 0.023 -0.089 -0.005 -0.039 -0.074 0.130 -0.105 -0.014 
0.004 0.081 0.056 -0.041 -0.099 


MEAN=-.267866E-02 VARIANCE = 0.198143E-01 


ENTER TITLE FOR PLOTS. 


autos and pautos of series c data / two nonseasonal differences 


YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 


Ley A 


0.144 
0.040 


0.117 


0.122 
-0.020 


0.094 





DO YOU WANT TO GO AGAIN? 
Y 


x 


ENTER LETTER FOR OPTION YOU WANT. 


C 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE:, &GOTO -QUES 

alter file ft02f£001 pl sere £t02f001 pl 
alter file £t03£001 pl file £t02£001 pl 
&goto -ques 


DO YOU WANT TO GO AGAIN? 
4 


ENTER LETTER FOR OPTION YOU WANT. 
e 
IS YOUR DATA IN FILE FTO2F001? 


y 
EXECUTION BEGINS... 
ENTER NUMBER OF AR PARAMETERS (EXCLUDE DIFFERENCES). 
0 
ENTER NUMBER OF MA PARAMETERS. 
2 
LENGTH OF TIME SERIES 2 224 
0 AR PARAMETERS 2 MA PARAMETERS 


MA PARAMETERS ARE: 
THETA(1)= 0.1378 


THETA(2)= 021296 

MA CONSTANT= -0.26787E-02 WHITE NOISE VARIANCE= Oel6I559e—0r 

CHI SQUARE RESIDUAL LACK OF FIT VALUE = 36.75 DF = 23 
SIGNIFICANCE = 0.0346 

DO YOU WANT TO PLOT AUTO AND PAUTO OF RESIDUALS? 


4 
*EXECUTION BEGINS... 
AUTOCORRELATIONS 


0.019 0.034 -0.134 -0.093 -0.012 -0.043 0.005 -0.063 -0.125 0.088 
-0.142 0.069 -0.098 0.025 0.004 -0.055 0.157 -0.090 0.005 -0.006 
0.009 0.133 0.005 -0.015 -0.133 


PARTIAL AUTOCORRELATIONS 
0.019 0.033 -0.135 -0.090 0.001 -0.055 =-0.019 -0.071 -0.143 0.087 
-0.166 0.017 -0.104 -0.013 =-0.025 -0.089 0.123 -0.123 -0.021 -0.025 
0.024 0.069 0.024 -0.069 =-0.102 
MEAN20.236684E=-02 VARIANCE = 0.189503E=-01 

ENTER TITLE FOR PLOTS. 

autos and pautos of residuals / arima(0,2,2) model of series c data 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 

PICK UP IN ROOM I140 UNDER YOUR USER ID NUMBER. 


DO YOU WANT TO GO AGAIN? 
Z 


ZS 





ENTER LETTER FOR OPTION YOU WANT. 


Cc 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE: &GOTO -QUES 

erase file £t02f001 

stat 

P (191): 27 FILES; 265 REC IN USE, 31 LEFT (of 296), 90% FULL (2 CYL) 
erase file ftO8f001 

alter sere £t02f001 pl file ft02f£001 pl 
&goto ~ques 

DO YOU WANT TO GO AGAIN? 

y 


ENTER LETTER FOR OPTION YOU WANT. 
d 
*ITS YOUR DATA IN FILE FTO2F001? 


y 
EXECUTION BEGINS... 
TS YOUR TIME SERIES SEASONAL? 
2 : 
ENTER NUMBER OF NONSEASONAL DIFFERENCES. 


a 
DO YOU WANT TO PLOT AUTO AND PAUTO OF TRANSFORMED DATA? 


n 
DO YOU WANT TO GO AGAIN? 
Y 


ENTER LETTER FOR OPTION YOU WANT. 


a4 

*ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE: &GOTO -QUES 

alter file £t02f001 pl sere ft02f001 pl 

alter file f£t03f001 pl file f£t02f001 pl 

&goto -ques 

DO YOU WANT TO GO AGAIN? 

y 


ENTER LETTER FOR OPTION YOU WANT. 
e 


{[S YOUR DATA IN FILE FTO2F001? 


sé 
EXECUTION BEGINS... 

ENTER NUMBER OF AR PARAMETERS (EXCLUDE DIFFERENCES). 
il 


ENTER NUMBER OF MA PARAMETERS. 

0 

* LENGTH OF TIME SERIES = 225 

1 AR PARAMETERS 0 MA PARAMETERS 


AR PARAMETERS ARE: 
PHI(1)= 0.8073 


MA CONSTANT= -0.66789E-02 WHITE NOISE VARIANCE= 0.177633E-01 

CHI SQUARE RESIDUAL LACK OP FIT VALUE | = 28. 88 DF = 24 
SIGNIFICANCE = 0.2247 ~~ oat 

DO YOU WANT TO PLOT AUTO AND PAUTO OF RESIDUALS? 


v4 
EXECUTION BEGINS... 
AUTOCORRELATIONS 


0.016 0.014 -0.049 -0.005 0.061 0.024 0.078 -0.018 -0.087 0.130 


~0.093 0.084 -~0.055 0.043 0.007 -0.040 0.169 -0.079 ~0.001 -0.010 
0.000 0.129 0.014 ~0.013 ~0.127 
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PARTIAL AUTOCORRELATIONS 


0.016 0.014 -0.050 -0.004 0.063 0.019 0.075 -0.015 -0.087 0.147 
-0.104 0.070 -0.048 0.042 0.003 -0.030 0.158 -0.088 0.013 -0.025 


MEAN#=-, 252517E=-02 VARIANCE = 0.177569E-01 


ENTER TITLE FOR PLOTS. 

autos and pautos of residuals / arima(1,1,0) model of series c¢ data 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 

PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 


DO YOU WANT TO GO AGAIN? 
Y 


* 
ENTER LETTER FOR OPTION YOU WANT. 
Cc 


ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE: &GOTO -QUES 

erase file f£t02f001 

alter serc £t02£001 pl file f£t02£001 pl 
&goto -ques 

DO YOU WANT TO GO AGAIN? 

y 


ENTER LETTER FOR OPTION YOU WANT. 
£ 

IS YOUR DATA IN FILE FTOQ2F001? 

Y 

EXECUTION BEGINS... 


IS YOUR SERIES SEASONAL? 
n 


ENTER NONSEASONAL DATA AS REQUESTED: 
ENTER NUMBER OF NON-SEASONAL DIFFERENCES. 

T . 

ENTER NUMBER OF NON-SEASONAL AR PARAMETERS. 
a. 


ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 
0 
NOW INPUT YOUR PARAMETER ESTIMATES: 


ENTER MA CONSTANT TERM, THETAO. 


ete oy-seasonat AR PARAMETER PHI(1). 

ee MAXIMUM FORECAST LEAD TIME VIA FORMAT I2. 
VENTER INDEX FOR PLOT ORIGIN VIA FORMAT 13. 

ENTER INDEX FOR FORECAST ORIGIN VIA FORMAT I3. 
“ENTER SIGNIFICANCE LEVEL FOR CONFIDENCE INTERVALS. 


WAS YOUR DATA TRANSFORMED IN THE TRANS PROGRAM? 
n : 


DO YOU WANT BASIC OUTPUT AT THE TERMINAL? 
Y 


IL Aas. 





THE LAST 10 WVEC VALUES: 


0.222000E 02 Q.218000E 02 O.213000E 02 0.208000E 02 
O.197000E 02 O.193000E 02 Q.191000E 02 O0.190000E 02 
THE 25 FORECAST VALUES: 
¢ 

0.203896E 02 Q.200517E 02 0O.197722E 02 0.195399E 02 
O.191821E 02 OQ.190434E O02 O0O.189248E 02 O.188223E 02 
O.186541E 02 Q.185837E 02 O0O.185203E 02 0.184623E 02 
O.183591E 02 OQ.I8312Z1E 02 O.182676E 02 Q.182249E 02 
O.181439E 02 Q.181050E 02 O.180669E 02 O0.180295E 02 
THE 25 UPPER FORECAST CONFIDENCE LIMITS: 

0.206094E 02 0.205056E 02 0.204779E 02 0.205033E 02 
O.206551E 02 O.207621E 02 O0.208812E 02 Q.210079E 02 
O.212721E 02 OQO.214056E 02 O.215381E 02 O.216687E 02 
O.219221E 02 O.220442E 02 O.221630E 02 Q.222783E 02 
0.224989E 02 0.226043E 02 0.227064E 02 Q.228055E 02 


ALPHA FOR THE CONFIDENCE LIMITS 1S: 0.100 


DO YOU WANT TO GO AGAIN? 


ME 
ENTER LETTER FOR OPTION YOU WANT. 


4 
EXECUTION BEGINS... 
SOLVES FOR ROOTS OF CHARACTERISTIC EQUATION. 


ENTER NO. OF AR PARAMETERS (UNDIFF. FORM). 


ENTER AUTOREGRESSIVE PARAMETERS: 
ENTER PHI(1). 


1.8073 
ENTER PHI(2). 
-0.8073 
THE 2 ROOTS OF CHARACTERISTIC EQ. AND NORMS ARE: 
eve 39 0.0 1.239 
1.000 0.0 1.000 


DO YOU WANT TO GO AGAIN? 
y 


ENTER LETTER FOR OPTION YOU WANT. 
ad 5 


*EXECUTION BEGINS... 
IS YOUR TIME SERIES SEASONAL? 
D 


ENTER NUMBER OF NONSEASONAL DIFFERENCES. 
0 . 


DO YOU WANT TO PLOT AUTO AND PAUTO OF TRANSFORMED DATA? 
n 


DO YOU WANT TO GO AGAIN? 
¥ 


ENTER LETTER FOR OPTION YOU WANT. 


Ss 
*EXECUTION BEGINS... 
ie NUMBER OF AR PARAMETERS (UNDIFFERENCED FORM) 


ENTER ESTIMATED AR PARAMETER PHI(1). 
oO 7S 


ENTER ESTIMATED AR PARAMETER PHI(2). 
-u.5073 


IL Aas 


0.202000E 
0.188000E 


0.193456E 
0.187329E 
0.184089E 
0.181838e 
0.179926E 


0.205662E 
0.211390E 
0.217969E 
0.223903E 
0.229016E 


02 
02 


02 
02 
02 
02 
02 


02 
02 
O02. 
02 
02 





AR PARAMETERS ARE 
ARE THESE OK? 


1.8073 


y é 
ENTER NUMBER OF MA PARAMETERS 


0 
ENTER OVERALL MA CONSTANT . 
-0.00667 


ao = 


‘ENTER ESTIMATED WHITE NOISE VAR 
0.0177 





-0.8073 


DO YOU WANT STARTING VALUES TO BE THE LAST VALUES 


OF THE ACTUAL SERIES? 


f 

START (1)= 19.000000 
START(2)= 18.799988 
STARTING’VALUE(S) OK? 


Y 
ENTER RANDOM NUMBER SEED (BETWEEN 0 AND 1). 


0.12345 


ENTER NUMBER OF VALUES YOU WANT TO SIMULATE VIA 13. 


50 


DO YOU WANT A PRINTOUT OF THE FIRST 20 VALUES? 


Y 
SIMULATED VALUES ARE: 


18.50224 18.33957 
17.85420 17.69398 
16.94655 16.68689 
16.31497 16.64629 

DO YOU WANT TO GO AGAIN? 

n 

CONTROL RSTURNED TO CMS 

R; ~ ras 


18.33694 
17.58267 
16.43227 


17.04015 


ey 


18.16435 
17.50722 
16.42499 
17..47137 


17.98804 
17.26422 
16.26595. 
17.62828 
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APPENDIX C 
KEYBOARD SESSION AND OUTPUT FROM THE ANALYSIS OF 
TIME SERIES G (SEASONAL) 

This appendix contains the complete terminal keyboard 
session and all offline output from the analysis of Bon 
and Jenkins’ [Ref. 4] time series G (International airline 
passengers, monthly totals, Jan 1949 - Dec 1960) using the 
Time Series Editor. This appendix is intended to supplement 
the discussion of the analysis of series G given in Chapter 
Mee Ot this report. 

The terminal listings are presented first, followed by 


the computer offline output. 
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repeat login nur@pn@ 


login 1621p44 400k 
ENTER PASSWORD: 


ENTER 4-DIGIT PROJECT NUMBER FOLLOWED BY 4-CHARACTER COST CENTER CODE: 
0444r172 

READY AT 15.12.34 ON 09/16/78 

CMS Version 3.25 


cP 
link 2069p 191 192 
ENTER PASSWORD: 


SET TO READ ONLY 


ipl cms 
CMS Version 3.25 


login 191 


R; 

login 192 t.p_ . ere es 
T (192) R/O : 

R; 


timeser 
EACH 2 SECONDS EXECUTION TIME IS INDICATED BY *. 


YOU HAVE ENTERED THE TIME SERIES SD ETON - 


» 4 ™ 


ee eS 


PLEASE RESPOND TO ‘BACH’ QUERY “WITH AN INPUT AT THE TERMINAL. 
ENTER ONLY THE FIRST LSTTER FOR A WORD RESPONSE. 

ENTER NUMERICAL VALUES VIA FORTRAN EEG 

TYPE INTEGER VALUES (RIGHT JUSTIFIED) FOR NAMES STARTING 

WITH I THRU N. TYPE FLOATING VALUES. WITH DECIMAL FOR ALL OTHERS. 


DO YOU WANT A LIST OF THE OPPIONS?. ov / 

4 * 

OPTION DESCRIPTION 

GENERATE ----- GENERATE ANY ARIMA TIMB SERIES © 

AUTO --------- CALCULATE AUTOCORRELATIONS, PAUTOS, MEAN AND VARIANCE 
PLOT -------=-PLOT A TIME SERIES 

ESTIMATE ----- CALCULATE MAX LIKELIHOOD ESTIMATES OF ARMA PARAMETERS 
DIFF --------- DIFFERENCE A TIME SERIES - 

FORECAST ----- FORECAST FUTURE VALUES, CONSTRUCT CONFIDENCE INTERVALS 
TRANS ------>s-TRANSFORMS VALUES OF A TIME SERIES 

ROOTS -------- DETERMINES ROOTS OF ARIMA CHARACTERISTIC EQUATION 
ZFORMAT ------ ALTER DATA FILE TO FORMAT 5P15 <6 

CMSWORK ------PERFORM CP/CMS COMMANDS IN TIMESER EXEC 

SIMULATE ----- SIMULATE NONSEASONAL TIME SERIES 

YESTSEAS =----CALCULATE INITIAL SEASONAL PARAMETERS 

WMARQRDT -----MARQUARDT SOLUTION FOR PARAMETER ESTIMATES 

XSUMSQ ------- CALCULATE SUM OF SQUARES FOR ARBITRARY PARAMETERS 


WOULD YOU LIKE MORE INFO?. _ : 3 a 


Yy ae AS wy 

ENTER OPTION YOU WANT INFO ABOUT. 

Y 

YESTSEAS ----THIS PROGRAM CALCULATES INITIAL PARAMETER ESTIMATES 


FOR SEASONAL AND NONSEASONAL ARIMA MODELS, TO BE USED AS INPUTS TO 
THE WMARQRDT PROGRAM. 


DO YOU WANT INFO ABOUT ANOTHER OPTION? 


y 
ENTER OPTION YOU WANT INFO ABOUT. 
a 


IL Sy; 





W4ARQRDT ----THIS PROGRAM CALCULATES NON=-LINEAR LEAST SQUARES 
ESTIMATES OF BOX-JENKINS PARAMETERS FOR A GENERAL (SEASONAL 
OR NONSEASONAL) ARIMA MODEL. IT REQUIRES INITIAL PARAMETER 
ESTIMATES AS STARTING VALUES; THESE MAY BE CALCULATED USING 
PROGRAM YESTSEAS. 


DO YOU WANT INFO ABOUT ANOTHER OPTION? 


Y 
ENTER OPTION YOU WANT INFO ABOUT. 
x 


fees ------ THIS PROGRAM CALCULATES THE SUM OF SQUARED RESIDUALS 
FOR ANY SET OF ARIMA (SEASONAL OR NONSEASONAL) PARAMETERS. 


DO YOU WANT INFO ABOUT ANOTHER OPTION? 


n 
*DO YOU WANT TO TRY A SESSION? 
y 


ENTER LETTER FOR OPTION YOU WANT. 


SG 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE: &GOTO -QUES 

stat 

P (191): 26 FILES; 235 REC IN USE, 61 LEFT (of 296), 79% FULL (2 CYL) 
Gog f 

FILES:s= NO RDR, NO PRT, NO PUN 

listf * £t02£001 

FILENAME FILETYPE MODE NO.REC. DATE 

SERG FTO2FOO] Pl 3 9/16 

SERC FTO2FOO] Pl 75 7 49716 

alter serg £t02f£001 pl file £t02f001 pl 
&goto -ques 

DO YOU WANT TO GO AGAIN? 


yi 2 
ENTER LETTER FOR OPTION YOU WANT. 


P 
IS YOUR DATA IN FILE FTO2F001? 


¥ 
EXECUTION BEGINS... 
ENTER TITLE FOR PLOT 
series g data / undifferenced 
TIME SERIES PLOTS HAVE BEEN PRINTED OFFLINE 
DO YOU WANT TO GO AGAIN? 
Y 


ENTER LETTER FOR OPTION YOU WANT. 
a 
IS YOUR DATA IN FILE FTO2F001? 


¥ 
EXECUTION BEGINS... 
AUTOCORRELATIONS 


0-948 0.876 0.807 0.753 0.714 0.682 0.663 0.656 0.671 0.703 
0.743 0.760 0.713 0.646 0.586 0.538 0.500 0.469 0.450 0.442 
0.457 0.482 0.517 0.532 0.494 

PARTIAL AUTOCORRELATIONS 
0.948 -0.229 0.038 0.094 0.074 0.008 0.126 0.090 0.232 0.166 


0.171 -0.135 -0.540 -0.027 0.091 0.025 0.032 0.0772 0.048 -0.046 
0.046 -0.100 0.052 0.048 -0.163 


Sal 





MEAN= 280.299 VARIANCE = 14292.0 


ENTER TITLE FOR PLOTS. 


autos and pautos of series g data / undifferenced 


YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 


DO YOU WANT TO GO AGAIN? 
Y 


E 
t 


NTER LETTER FOR OPTION YOU WANT. 


IS YOUR DATA IN FILE FTO2F001? 


y- 
EXECUTION BEGINS... 


DO YOU WANT-TO TRANSLATE THE ORIGIN: 


nN 


DO YOU WANT TO RESCALE THE VALUES: 


nN 


W=Z-SHIFT? 


W2Z *SCALE? 


DO YOU WANT A LOG TRANSFORMATION? : 


Y 
TRANSFORMATION IS W(I)#LOG(SCALE*(2Z(1I)-SHIFT)-FACTOR) WHERE: 


SCALE= 1.00000 SHIFT=0.0 FACTOR20.0 


YOUR ORIGINAL SERIES IS IN DATA FTO2F001 Pl. 
YOUR TRANSFORMED SERIES IS IN FILE FTO2F001 Pl. 


DO YOU WANT TO PLOT THE TRANSFORMED VALUES? 


Y 
EXECUTION BEGINS... 


* 


P 


ENTER TITLE FOR PLOT 
lot of series g data / natural log transform 


TIME SERIES PLOTS HAVE BEEN PRINTED OFFLING 


DO YOU WANT TO GO AGAIN? 


Y 


ENTER LETTER FOR OPTION YOU WANT. 


a 


IS YOUR DATA IN FILE FTO2F001? 


Y 
EXECUTION BEGINS... 


a 


D 
Y 


AUTOCORRELATIONS 
0.954 0.899 0.851 0.808 0.779 0.756 0.738 0.727 0.734 
0.758 0.762 0.717 0.663 0.618 0.576 0.544 0.519 0.501 
0.498 0.506 0.517 0.520 0.484 - 
PARTIAL AUTOCORRELATIONS 
0.954 -0.118 0.054 0.024 0.116 0.044 0.038 0.100 0.204 
0.106 -0.042 -0.485 -0.034 0.042 -0.044 0.028 0.037 0.042 
0.073 -0.033 0.061 0.031 ~0.194 . 
MEAN= 5.54218 VARIANCE = 0.193531 


ENTER TITLE FOR PLOTS. 

utos and pautos of series g data / natural log transform 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 

O YOU WANT TO GO AGAIN? 


15) 2 


0.744 
0.490 


0.064 
0.014 





ENTER LETTER FOR OPTION YOU WANT. 
| 
IS YOUR DATA IN FILE FTO2F001? 


of 
EXECUTION BEGINS... 
IS YOUR TIME SERIES SEASONAL? 


‘4 
ENTER ORDER OF SEASONAL DIFFERENCING. 


ik 
ENTER LENGTH OF SEASONAL PERIOD VIA I2. 
v2 
ENTER NUMBER OF NONSEASONAL DIFFERENCES. 


1 
DO YOU WANT TO PLOT AUTO AND PAUTO OF TRANSFORMED DATA? 


y 
EXECUTION BEGINS... 
AUTOCORRELATIONS 


“0.341 0.105 -0.202 0.021 0.056 0.031 -0.056 -0.001 0.176 -0.076 
0.039 -0.091 0.223 -0.018 -0.100 


PARTIAL AUTOCORRELATIONS 


-0.341 -0.013 -0.193 =r aS 0.033 0.035 -0.060 -0.020 0.226 0.043 
0.047 -0.339 -0.109 -0.077 -0.022 -0.140 0.026 0.115 -0.013 -0.167 
0.132 -0.072 0.143 -0.067 -0.103 ‘ , 


MEAN=0 .290920E-03 VARIANCE 2 0.208604E-02 


ENTER TITLE FOR PLOTS. 

autos and pautos for series g s@data / 1 ns diff, 1 seas diff, ln xform 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 

PICK UP IN ROOM I140 UNDER YOUR USER ID NUMBER. 


DO YOU WANT TO GO AGAIN? 
Y 


ENTER LETTER FOR OPTION YOU WANT. 

Cc 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE: &GOTO -QUES 

alter file £t03f001 pl save £t03f£001 pl 
alter data ftO2f001 pl serg ft02f001 pl 

Stat 

P (191): 29 FILES; 267 REC IN USE, 29 LEFT (of 296), 90% FULL (2 CYL) 
erase file ftO08f00l 

ao g f 

&goto -ques 

DO YOU WANT TO GO AGAIN? 

Y 


ENTER LETTER FOR OPTION YOU WANT. 

ts YOUR DATA IN FILE FT02F001? 

EXECUTION BEGINS... 

erer LENGTH OF SEASON VIA FORMAT 12. 
etek NUMBER OF NON-SEASONAL DIFFERENCES. 
mrer NUMBER OF SEASONAL DIFFERENCES. 


_ °* NUMBER OF NON-SEASONAL AR PARAMETERS. 


oS 





ENTER NUMBER OF SEASONAL AR PARAMETERS. 


0 
ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 


1 : 
ENTER NUMBER OF SEASONAL MA PARAMETERS. 


1 


INITIAL PARAMETER ESTIMATES FOR MARQRT : 


THETA(1) = 
THETAS (1) # 


0.390425 


0.534397 


DO YOU WANT TO GO AGAIN? 


Y 


ENTER LETTER FOR OPTION YOU WANT. 


W 


IS YOUR DATA IN FILE FTO2F001? 


Y 
EXECUTION BEGINS... 


ENTER 


1 
ENTER 
1 
ENTER 
0 
ENTER 
0 
ENTER 
1 
ENTER 
1 
ENTER 
12 


NUMBER 
NUMBER 
; 
NUMBER 
NUMBER 
NUMBER 
NUMBER 


LENGTH 


OF 


OF 


OF 


OF 


OF 


OF 


OF 


NON-SEASONAL DIFFERENCES. 
SEASONAL DIPFERENCES. 
NON-SEASONAL AR PARAMETERS. 
SEASONAL AR PARAMETERS. 
NON-SEASONAL MA PARAMETERS. 
SEASONAL MA PARAMETERS. 


SEASON VIA FORMAT I2. 


NOW INPUT YOUR INITIAL PARAMETER ESTIMATES, AS REQUESTED. 


ENTER NON-SEASONAL MA PARAMETER THETA(1). 


~390425 


ENTER SEASONAL MA PARAMETER THETAS(1). 


~ 934397 


x 


CONVERGENCE HAS BEEN REACHED IN MAX( 7, 


7) .(QITERATIONS. 


SELECTED OUTPUT FOLLOWS: 


PARAMETER ESTIMATE 


THETA(1) = 
THETAS(1) = 


MOVING AVERAGE CONSTANT: THETAO 2 


STANDARD ERROR 


06377152 
G.572387 


0.819933D-01 
0.780252D-01 


0.000291 


CHI-SQUARE STATISTIC FOR RESIDUAL LACK OF FIT 2 29.713571 


DEGREES OF FREEDOM = 38 


PROBABILITY OF EXCEEDING STATISTIC = 


0.829438 


154 


DO YOU WANT TO PLOT AUTO AND PAUTO OF RESIDUALS? 
y 


EXECUTION BEGINS... 
AUTOCORRELATIONS 


0.009 0.026 -0.129 -0.105 0.078 0.077 -0.036 -0.033 0.103 -0.050 
0.026 -0.020 0.013 0.035 0.067 -0.130 0.052 0.015 -0.093 -0.090 
-0.027 -0.014 0.213 0.009 -0.042 


PARTIAL AUTOCORRELATIONS 


> 


0.009 0.026 -0.130 -0.104 0.088 0.068 -0.072 -0.029 0.151 -0.057 
-0.022 0.018 0.046 0.001 0.053 -0.114 0.069 0.029 -0.129 -0.134 
0.042 -0.014 0.143 -0.025 0.013 


MEAN=0.236005E-02 VARIANCE = 0.139713E-02 


ENTER TITLE FOR PLOTS. 

autos and pautos of residuals / (0,1,1)x(0,1,1)12 model of series g 
YOUR AUTO AND PAUTO PLOTS HAVE BEEN PRINTED OFFLINE. 

PICK UP IN ROOM 1140 UNDER YOUR USER ID NUMBER. 

® . 


DO YOU WANT TO GO AGAIN? 
y E 


ENTER LETTER FOR OPTION YOU WANT. 
4 ¥ e 
IS YOUR DATA IN FILE FTO2F001? 


y 

EXECUTION BEGINS... 
ENTER LENGTH OF SEASON VIA FORMAT [2. 

12 : 
ENTER NUMBER OF NON-SEASONAL DIFFERENCES. 

l 


ENTER NUMBER OF SEASONAL DIFFERENCES. 
Perce NUMBER OF NON=-SEASONAL AR PARAMETERS. 
Pirer NUMBER OF SEASONAL AR PARAMETERS. 
Pree NUMBER OF NON-SEASONAL MA PARAMETERS. 
“ENTER NUMBER OF SEASONAL MA PARAMETERS. 


rk 


INITIAL PARAMETER ESTIMATES FOR MARQRT 3: 


PHI(1) = 0.112699 
THETA(1) = 0.490883 
THETAS(1) = 0.533679 


DO YOU WANT TO GO AGAIN? 
Y 


ENTER LETTER FOR OPTION YOU WANT. 
W 
IS YOUR DATA IN FILE FT02F001? 


af 
EXECUTION BEGINS... 
_* NUMBER OF NON-SEASONAL DIFFERENCES. 


ENTER NUMBER OF SEASONAL DIFFERENCES. 
1 


ENTER NUMBER OF NON=SEASONAL AR PARAMETERS. 
1 


Rope, 





ENTER NUMBER OF SEASONAL AR PARAMETERS. 


0 

ENTER NUMBER OF NON=-SEASONAL MA PARAMETERS. 
1 

ENTER NUMBER OF SEASONAL MA PARAMETERS. 


1 
ENTER LENGTH OF SEASON VIA FORMAT I2. 
7 


NOW INPUT YOUR INITIAL PARAMETER ESTIMATES, AS REQUESTED. 


ENTER NON-SEASONAL AR PARAMETER PHI(1). 
0.112699 

ENTER NON~SEASONAL MA PARAMETER THETA(1). 
0.490883 

ENTER SEASONAL MA PARAMETER THETAS(1). 
0.533679 


eek 


CONVERGENCE HAS BEEN REACHED IN MAX(10,10) ITERATIONS. 


SELECTED OUTPUT FOLLOWS: 


PARAMETER ESTIMATE STANDARD ERROR 
PHT aad 146667 0.224464D 00 

- ° THETA(1) = 0.509816 0.195934D 00 
THETAS(1) = 0.573388 0.784896D-01 


MOVING AVERAGE CONSTANT: THETAO = 0.000248 


CHI-SQUARE STATISTIC FOR RESIDUAL LACK OF FIT = 28.302543 
DEGREES OF FREEDOM = 37 


PROBABILITY OF EXCEEDING STATISTIC = 0.847102 
DO YOU WANT TO PLOT AUTO AND PAUTO OF RESIDUALS? 


y 
*EXECUTION BEGINS... 
AUTOCORRELATIONS 


“0.013 0.073 -0.100 -0.088 0.077 0.073 -0.031 -0.033 0.105 -0.054 
0.033 -0.018 0.014 0.027 0.070 -0.133 0.051 0.008 -0.092 -0.084 
-0.020 -0.023 0.208 0.005 -0.036 


? 


PARTIAL AUTOCORRELATIONS 


“0.013 0.073 -0.098 -0.096 0.092 0.081 -0.064 -0.040 0.152 -0.050 
“0.027 0.017 0.047 -0.004 0.053 -0.118 0.055 0.031 -0.123 -0.135 
0.051 -0.006 0.148 -0.019 0.003 


MEAN=0.253052E-02 VARIANCE 2 0.139373E-02 
ENTER TITLE FOR PLOTS. : 
autos and pautos of residuals / (1,1,1)x(0,1,1)12 model of series g 
YOUR AUTO AND PAUTO PLOTS'‘HAVE BEEN PRINTED OFFLINE. 
PICK UP IN ROOM I140 UNDER YOUR USER ID NUMBER. 


DO YOU WANT TO GO AGAIN? 
y 


IL SYS 





mers 


ENTER LETTER FOR OPTION YOU WANT. 


£ 


IS YOUR DATA IN FILE FTO2F001? 


Y 
EXECUTION BEGINS... 


IS YOUR SERIES SEASONAL? 


Y 


ENTER SEASONAL DATA AS REQUESTED: 


- ENTER 
i 
ENTER 
0 
ENTER 
1 


ENTER 
12 


ENTER 
_ ENTER 
1 

ENTER 


0 
ENTER 
it 


NUMBER OF SEASONAL DIFFERENCES. 


LENGTH OF SEASON VIA FORMAT 12. 


NONSEASONAL DATA AS REQUESTED: 


NOW INPUT YOUR PARAMETER ESTIMATES: 


ENTER MA CONSTANT TERM, THETAO. 


~000291 


ENTER NON~SEASONAL MA PARAMETER THETA(1). 


ee 7 152 


NUMSER OF SEASONAL AR PARAMETERS. 


NUMBER OF SEASONAL MA PARAMETERS. 


NUMBER OF NON~SEASONAL DIFFERENCES. 


ENTER SEASONAL MA PARAMETER THETAS(1). 


372387 


NUMBER OF NON~SEASONAL AR PARAMETERS. 


NUMBER OF NON-SEASONAL MA PARAMETERS. 


ENTER MAXIMUM FORECAST LEAD TIME VIA FORMAT I2. 


30 


ENTER INDEX FOR PLOT ORIGIN VIA FORMAT I3. 


100 


ENTER INDEX FOR FORECAST ORIGIN VIA FORMAT I3. 


rl 


ENTER SIGNIFICANCE LEVEL FOR CONFIDENCE INTERVALS. 


0.10 


WAS YOUR DATA TRANSFORMED IN THE TRANS PROGRAM? 


4 


DO YOU WANT BASIC OUTPUT AT THE TERMINAL? 


“4 


THE LAST 10 WVEC VALUES: 


0.418999E 
0.606000E 


03 
03 


0.461001E 


03 


0.507999E 03 


THE 30 FORECAST VALUES: 


0.398539E 
0.469710E 
0.460421E 
0.520593E 
0.705666E 
0.526812E 


03 
03 
03 
03 
03 
03 


O.415184E 
0.543456E 
0.405213E 
0.507548E 
0.590927E 
0.502002E 


03 
03 
03 
03 
03 
03 


0.472000E 
0.461001E 


0.395402E 
0.618194E 
0.448016E 
G25267 715 
0.519089E 
0.587782E 


4b 5) 7 


03 
03 


03 
03 
03 
03 
03 
03 


0.535001E 
0.390001& 


0.462697E 
0.626276E 
0.466864E 
O.611991E 
0.456978E 
0.573220E 


03 
03 


03 
03 
03 
03 
03 
03 


0.622000E 
0.432002E 


0.450971E 
0.524292E 
0.444748E 
0.696357E 
0.505397E 
Ue 2973005 


03 
03 


03 
03 
a 
03 
03 
03 





THE 30 UPPER 


0.399598E 
0.470813E 
0.461557E 
O.521779E 
0.706894E 
0.528090E 


03 
03 
03 
03 
03 
03 


0.416254E 
0.544567E 
0.406354E 
0.508743E 
0.592163£ 
0.503292E 


03 
03 
03 
03 
03 
03 


0.396481E 
0.619311E 
0.449173E 
0.529995E 
0.520332E 
0.589083E 


FORECAST CONFIDENCE LIMITS: 


03 
03 
03 
03 
03 
03 


0.463785E 
0.627399E 
0.468030£E 
0.613203E 
0.458229E 
0.574532E 


ALPHA FOR THE CONFIDENCE LIMITS IS: 0.100 


* 


DO YOU WANT TO GO AGAIN? 


Y 


ENTER LETTER FOR OPTION YOU WANT. 


3 


Is YOUR DATA IN FILE FTO2F001? 


uf 
EXECUTION BEGINS... 


IS YOUR SERIES SEASONAL? 


Y 
ENTER 


ENTER 
i: 

ENTER 
0 

ENTER 
1 

ENTER 
12 


ENTER 
ENTER 


1 
ENTER 
1 


SEASONAL DATA AS REQUESTED: 


NUMBER OF SEASONAL 


¢ 


NUMBER OF SEASONAL AR PARAMETERS. 


DIFFERENCES. 
NUMBER OF SEASONAL MA PARAMETERS. 


LENGTH OF SEASON VIA FORMAT I2. 


NONSEASONAL DATA AS REQUESTED: 
NUMBER OF NON=SEASONAL DIFFERENCES. 


NUMBER OF NON=SEASONAL AR PARAMETERS. 


ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 


NOW INPUT YOUR PARAMETER ESTIMATES: 


ENTER MA CONSTANT TERM, THETAO. 
.000248 


ENTER NON-SEASONAL AR PARAMETER PHI(1). 
0.146667 

ENTER NON=SEASONAL MA PARAMETER THETA(1). 
0.509816 


ENTER SEASONAL MA PARAMETER THETAS(1). 
0.573388 


ENTER MAXIMOM FORECAST LEAD TIME VIA FORMAT I2. 
30 


ENTER INDEX FOR PLOT ORIGIN VIA FORMAT I3. 
100 


ENTER INDEX POR FORECAST ORIGIN VIA FORMAT I3. 
131 


ENTER SIGNIFICANCE LEVEL FOR CONFIDENCE INTERVALS. 
0.10 


WAS YOUR DATA TRANSFORMED IN THE TRANS PROGRAM? 
Y 


DO YOU WANT.BASIC OUTPUT AT THE TERMINAL? 
“f 


leases: 


03 


03 
03 
03 
03 


0.452068E 
0.525422E 
0.445925E 
O0.697577E 
0.506664E 
0.598708E 


03 
03 
03 
03 
03 
Q2 





THE LAST 10 WVEC VALUES: 


0.418999E 
0.606000E 


0.398650, 


0.469348E 
0.460071E 


(0.520096E 


0.704905E 
0.526113E 


03 
03 


03 
03 
03 
03 
03 
03 


THE 30 UPPER 


0.461001E 
0.507999E 


. THE 30 FORECAST VALUES: 


0.414875E 
0.543045E 
0.404899E 
0.507066E 
0.590336E 
0.501400E 


03 
03 


03 
03 
03 
03 
03 
03 


0.472000E 
0.461001E 


0.395125E 
0.617700E 
0.447739E 
0.528257E 
0.518566E 
0.587074E 


IS YOUR DATA IN FILE FTO2F001? 


y 
EXECUTION BEGINS... 


IS YOUR SERIES SEASONAL? 


y 
ENTER SEASONAL DATA AS REQUESTED: 


1 
0 


ENTER 
a4 


1 


1 
0 
1 


Sd 


ENTER 


ENTER 


ENTER 
2 


ENTER 
ENTER 


ENTER 


NUMBER OF SEASONAL DIFFERENCES. 


LENGTH OF SEASON VIA FORMAT I2. 


NONSEASONAL DATA AS REQUESTED: 


03 
03 


03 
03 
03 
03 
03 
03 


FORECAST CONFIDENCE LIMITS: 


03 
03 
03 
03 
03 
03 


0.399709E 03 0.415945E 03 0.396204E 
0.470448E 03 0.544151E 03 0.618812E 
0.461199E 03 0.406032E 03 0.448887E 
0.521270E 03 0.508249E 03 0.529448E 
0.706117E 03 0.591556E 03 0.519792E 
0.527371E 03 0.502668E 03 0.588353E 
ALPHA FOR THE CONFIDENCE LIMITS IS: 0.100 
DO YOU WANT TO GO AGAIN? 

y 

ENTER LETTER FOR OPTION YOU WANT. 

£ 


NUMBER OF SEASONAL AR PARAMETERS. 


NUMBER OF SEASONAL MA PARAMETERS. 


NUMBER OF NON-SEASONAL DIFFERENCES. 


0.535001E 
0.390001E 


0.462366E 
0.625753E 
0.466360E 
0.611381E 
0.456512E 
O2.572532E 


0.463452E 
0.626870E 
0.467517E 
0.612579E 
0.457745E 
0.573821E 


NUMBER OF NON=-SEASONAL AR PARAMETERS. 


ENTER NUMBER OF NON=SEASONAL MA PARAMETERS. 


NOW INPUT YOUR PARAMETER ESTIMATES: 


ENTER MA CONSTANT TERM, 


000291 


ENTER NON-SEASONAL MA PARAMETER THETA(1). 


377152 


ENTER SEASONAL MA PARAMETER THETAS(1). 


972387 


THETAO. 


ENTER MAXIMUM FORECAST LEAD TIME VIA FORMAT I2. 
40 


L539 


03 
03 


03 
03 
03 
03 
03 
03 


03 
03 
03 
03 
03 
03 


0.622000E 
0.432002E 


0.450651E 
0.523897E 
0.444325E 
0.695632E 
0.504960E 
0.596632E 


0.451744E 
0.525020E 
0.445491E 
0.696837E 
0.506207E 
0.597930E 


03 
03 


03 
03 
02 
03 
03 
03 


03 
03 
03 
03 
03 
03 


ENTER INDEX FOR PLOT ORIGIN VIA FORMAT I3. 

075 

ENTER INDEX FOR FORECAST ORIGIN VIA FORMAT 13. 

120 

ENTER SIGNIFICANCE LEVEL FOR CONFIDENCE INTERVALS. 
QO. 10 


WAS YOUR DATA TRANSFORMED IN THE TRANS PROGRAM? 
sf 


DO YOU WANT BASIC OUTPUT AT THE-TERMINAL? 


YOU WANT TO GO AGAIN? 


eS as 


ENTER LETTER FOR OPTION YOU WANT. 
d 
IS YOUR DATA IN FILE FTO2F001? 


y 
EXECUTION BEGINS... 
IS YOUR TIME SERIES SEASONAL? 


, ENTER ORDER OF SEASONAL DIFFERENCING. 
1 


ENTER LENGTH OF SEASONAL PERIOD VIA 12. 
12 

ENTER NUMBER OF NONSEASONAL DIFFERENCES. 
a4 


DO YOU WANT TO PLOT AUTO AND PAUTO OF TRANSFORMED DATA? 


n 
DO YOU WANT TO GO AGAIN? 
Y 


ENTER LETTER FOR OPTION YOU WANT. 
2 ; 

ENTER DESIRED CP/CMS COMMANDS, ONE PER LINE. 
WHEN FINISHED TYPE: &GOTO -QUES 

alter file £t02£001 pl lnserg ft02f001 pl 
alter file £t03f001 pl file ft0O2f001 pl 
&goto -ques 

DO YOU WANT TO GO AGAIN? 

Y 


ENTER LETTER FOR OPTION YOU WANT. 
x 
IS YOUR DATA IN FILE FTO2F001? 


4 
EXECUTION BEGINS... 


IS YOUR SERIES SEASONAL? 
“ENTER LENGTH OF SEASON VIA FORMAT I2. 
Sener NUMBER OF SEASONAL AR PARAMETERS. 
ENTER NUMBER OF SEASONAL MA PARAMETERS. 
ENTER NUMBER OF NON=-SEASONAL AR PARAMETERS. 


ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 


~ Oo fF © 


NOW INPUT YOUR INITIAL PARAMETER ESTIMATES, AS REQUESTED. 


ENTER NON-SEASONAL MA PARAMETER THETA(1). 
0377152 


ENTER SEASONAL MA PARAMETER THETAS(1). 
~ 572387 


160 





FOR THE GIVEN INPUT PARAMETERS AND MODEL 


THE RESIOUAL SUM OF SQUARES IS: 0.18192935D 00 


DO YOU WANT TO TEST OIFFERENT PARAMETER VALUES? 
n 


DO YOU WANT TO TEST A DIFFERENT MODEL? 
Y 
IS YOUR SERIES SEASONAL? 
Y 
ENTER LENGTH OF SEASON VIA FORMAT I2. 
12 
ENTER NUMBER OF SEASONAL AR PARAMETERS. 
0 
ENTER NUMBER OF SEASONAL MA PARAMETERS. 
it 
ENTER NUMBER OF NON-SEASONAL AR PARAMETERS. 
1 
ENTER NUMBER OF NON-SEASONAL MA PARAMETERS. 
al 
NOW INPOT YOUR INITIAL PARAMETER ESTIMATES, AS REQUESTED. 


ENTER NON-SEASONAL AR PARAMETER PHI(1). 


- 146667 

ENTER NON-SEASONAL MA PARAMETER: THETA(1). 
- 509816 

ENTER SEASONAL MA PARAMETER THETAS(1). 
. 9/3388 


t 


FOR THE GIVEN INPUT PARAMETERS AND MODEL 


THE RESIDUAL SUM OF SQUARES IS: 0.18164378D 00 


DO YOU WANT TO TEST OIFFERENT PARAMETER VALUES? 
n 


DO YOU WANT TO TEST A DIFFERENT MODEL? 
n 


DO YOU WANT TO GO AGAIN? 
n 

CONTROL RETURNED TO CMS 
Rs; ; 


How 


QIINSAYSASTONN 7 VYivI 9 SSIUIS 
SLINA TO 5SE98°O =e SS IVIS “A 
SLINAN TO S6LT°O =wre s3¥V9OS-xX 
0I° oT Gc°sotT O0S°2el GL°9E CO °T 


REFEREERERERE EERE ERE REE EREREERFERREREE RE EERE EERE RERREREE HERERO RHE EEE TREE RH EREE Ot 


+ + 
“ e @ e * 
os t e @ ee & 
* e e 6a f& 
e e ee 
od e oe @ ee : 
> x ; « 
& e e td] « 
« e 6 ® 6 € 
* e ee e ee * 
+ a e @e e $ 
? 4" { ® e ® 6 « 
* ; nm e @ rs & 
* | : # 
rd | ‘ ee 6 e & 
ax @ e @e ® ee e * 
« ry e ® ® ® 
* ° * 
* - ¢ 
* - e eee a ee . 
+ e ee 4 
* ae be * 
* ' bd « 
1 { ° ooo i) € 
me e e H e ee ¢€ 
* ; e | ee * 
« * 
* be . * 
ny * e e ® e : * 
* ee ® a ’ 
+ @ e @ @@ e + 
* : x 
* : ’ 
* s % « 
+ ° nS 
* ee . & e » A ¢ 
* ee e e % 
« . es ‘ ' 
; : * 
* * 
“ . 
coe |e ¢ 
@ @@ P" 
* a ° + 
* * 
€ : ke 
* . * 
* pd * 
* ¢€ 
* * 
+ © + 
* : ” 
x * 
* : k 
* \ a 
* * 
* « 
* « 
*x . * 
* * 
+ © + 
KOE EOD ORE FERRER ERE KE HERE ERE RAKE EEE he eR RRS RR KEKE REE KES EEK K RE OE KR 
20 °ooT $2°80T OS °cl GL°9E 03° I 


00 °4DT 


ce°OStT 


L9°9L2 


00°: 92 


CE°S 44 


L9°SES 


22°29 


GZ 


ALTCCORRELAT ICNS 


TUWustw 
Oorar 
ror Alet 
e® ep &© @ 


Coccoe 


omd 7) LN et 
FHL 2d ot 
OTM 
e@e @ 8 


CwWOO 


QF Dep 
nO fas 
OPT May 
oeee 


aces 


MODM 
oh oh @) 
OUrene 
e®%® @e 8 


oocoe 


NOOO 
Oman 
OWS 
ee 8 @ 
OWwWo 


Pros 
at OW OY 
MAST 
@et%e 


COoOoO 


A ywiwge 
Ve tod he be 
fre OLN 
eo* @¢ @ 


Oo0e 


tn | ed dae! 
Onn 
COP LAY 
eo @ @ e 


“voco% 


WON? 
m0 OO 
hd fee Pls 
eee @ 


WOWOD 


Dey rr wad 
Fie ult 
Ape TN 
eoeoeoe 


Oooooe 


PARTIAL AUTOCORRELATIONS 


VON 
Of Himd 


NI 00.0 
VO oa 
NOVO 
ee @ @ 
fol olele) 


omvon 
aH OO 
Qeroo 
@®es @ 


eee 


ONO 
moos 


163 


14292.0 


VARIANCE 2 
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APPENDIX D 


PROGRAM LISTINGS 


This Appendix contains program listings of all programs 
resident in the Time Series Editor. It does not include 
listings of subroutines used from other libraries, such 


as IMSL routines. 
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APPENDIX E 
LISTING OF DATA SETS RESIDING 
[TN LHe, LIME SEREES EBITOR 
This appendix contains listings of the data sets now 

residing in the Time Series Editor. The data consists of 
Box and peeing ' [Ref. 4] time series C [Ref. 4, p. 528] 
and time series G [Ref. 4, p. 531], as well as a data set 
containing monthly Monterey, California, rainfall data 


from January 1931 through December 1976 . 
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